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This review article discusses the experimental 
and theoretical status of partonic charge 
symmetry. It is shown how the partonic content 
of various structure functions gets redefined when 
the assumption of charge symmetry is relaxed. 
We review various theoretical and 
phcnomenological models for charge symmetry 
violation in parton distribution functions. We 
summarize the current experimental upper limits 
on charge symmetry violation in parton 
distributions. A series of experiments are 
presented, which might reveal partonic charge 
symmetry violation, or alternatively might lower 
the current upper limits on parton charge 
symmetry violation. 



The notion of isotopic spin was introduced to 
account for the strong similarity between the pro- 
ton and neutron. In this picture the proton and 
neutron are defined as two components of a sin- 
gle object, the nucleon. If the strong interaction 
H s does not distinguish between the proton and 
neutron; then H s will commute with the isospin 
vector T, such that 

[H S ,T} = 0. (1) 

A strong interaction satisfying Eq. ([IJ is said to 
satisfy charge independence. Charge symmetry is 
a specific operation involving the isospin vector. 
It is defined as a rotation of 180° about the "2" 
axis in isospin space. Thus the charge symmetry 
operator P cs is defined as 

Pes = . (2) 

Charge symmetry involves interchanging a proton 
and neutron. When operating on light quarks, 
the charge symmetry operator interchanges up and 
down quarks, namely 

P os \u) = -\d) ; P cs \d) = \u). (3) 

In nucleon isospin space, the operation of charge 
symmetry thus interchanges up and down quarks 
(also up and down antiquarks), while interchang- 
ing proton and neutron labels. 

As is well understood, inclusive processes at 
high energies can be described in terms of a small 
number of structure functions, and these structure 
functions can be characterized in terms of parton 
distribution functions, or PDFs, that describe the 
probability of finding a given flavor quark or anti- 
quark with a fraction x of the nucleon 's momen- 
tum. Over the past 30 years, increasingly precise 
measurements have been made of the PDFs and 
their dependence on x and Q 2 . If we assume that 
charge symmetry is obeyed at the level of parton 
distributions, this implies the relations 

u p {x,Q 2 ) = d'\x,Q 2 ); 

d p (x,Q 2 ) = u n (x,Q 2 ); 

s p (x,Q 2 ) = s n (x,Q 2 ) = s(x,Q 2 ); 

c p (x,Q 2 ) = c n (x,Q 2 ) = c(x,Q 2 ). (4) 

In Eq. Q , the superscript describes the target nu- 
cleon and the quantities u, d, s and c represent the 



flavor of the struck quark. Relations analogous to 
Eq. Q are obtained by replacing all quark distri- 
butions by antiquarks. 

Until recently, all quark/parton phcnomcnolog- 
ical models assumed the validity of charge sym- 
metry at the outset. This was a sensible assump- 
tion for several reasons. First, charge symmetry 
is obeyed to a very high precision at low energies; 
whereas in many nuclear reactions isospin symme- 
try is obeyed only to the level of a few percent, in 
most cases charge symmetry is valid to better than 
one percent [JJ [2], Recent precise measurements 
of charge symmetry violation (CSV) in single-pion 
production in few-body systems [21 H] have led to 
better understanding of charge symmetry at low 
energies El EJ [7J EJ E] • The high precision of charge 
symmetry at low energies makes it natural to as- 
sume that charge symmetry is valid at high en- 
ergies; indeed, it is difficult to imagine a scenario 
with large charge symmetry violation at the par- 
tonic level, which would lead to very small CSV 
at low energies. Second, the assumption of charge 
symmetry reduces by a factor of two the number 
of independent quark PDFs that must be deter- 
mined. Third, early measurements of high en- 
ergy structure functions showed that the require- 
ments of charge symmetry were at least qualita- 
tively obeyed [TO] [TTJ [TJJ [131 El Ell- 
in 1998 the current situation regarding parton 
charge symmetry was reviewed |16j . Since then 
there have been several developments that war- 
rant an updated review. At that time, compar- 
ison of the F% structure functions from charged 
lepton deep inelastic scattering (DIS) and neu- 
trino charge-changing DIS [T7l [18] suggested sub- 
stantial CSV contributions in the nuclcon sea 
[Ti?l |2"U] , However, re-analysis of the neutrino re- 
actions [HI [22] removed the discrepancies that ap- 
peared at that time to indicate the possibility of 
surprisingly large CSV effects [19] . 

In the past few years CSV terms have for the 
first time been included in global fits to high en- 
ergy data [23 . Although these global fits contain 
some model dependence, nonetheless such fits al- 
low one to set phenomenological limits on CSV 
contributions to PDFs. In addition, another mech- 
anism for isospin violation in PDFs (quantum elec- 
tromagnetic or "QED splitting" effects) has now 
been included in calculations of PDFs (24] [25] 
through modification of what is termed DGLAP 
evolution [2U H7J HH]- Inclusion of these QED 
splitting terms also leads to CSV effects in par- 
ton distribution functions. 

The phenomenological limits for PDFs obtained 



from global fits to high energy data provide ef- 
fective upper limits for the magnitude of CSV ef- 
fects. As we shall see, these limits are somewhat 
larger than those obtained from theoretical esti- 
mates of partonic CSV contributions. Using these 
phenomenological estimates provides limits to the 
size of CSV effects that might reasonably be ob- 
served in certain experiments. We will use the 
phenomenological limits obtained from global fits 
to estimate the maximum value of CSV effects that 
might be seen in dedicated experiments. This will 
provide at least qualitative estimates of the size 
of effects that could be observed in various experi- 
ments. It will also provide guidance as to the most 
promising experiments that could tighten the ex- 
isting upper limits on parton CSV. 

Since the publication of our previous review 
on parton CSV, the NuTeV group has measured 
total cross sections for v and V charged-current 
and neutral-current reactions on an iron target 
[29l [30] , These measurements allow them to ex- 
tract an independent measurement of the Wein- 
berg angle. Their measurement differs by 3er from 
the extremely precise values for the Weinberg an- 
gle measured at the Z° mass [31] . The publication 
of the NuTeV measurement has led to a great deal 
of investigation of effects that might explain this 
result. We will review various 'QCD corrections' 
to the NuTeV result (i.e., corrections within the 
Standard Model) and in particular we will summa- 
rize the potential corrections to the NuTeV mea- 
surement from parton CSV. 

Our review is organized as follows. In Sect. |II A| 
we review the general form of high-energy cross 
sections in terms of structure functions. In Sect. 



II B we define parton distributions when charge 
symmetry violation is included. In Sect. |IIC| 
we review the definitions of structure functions in 
terms of quark/parton distributions. We list the 
general form of structure functions when one re- 
laxes the assumption of charge symmetry. In Sect. 
|IID| we write down relations between leading-order 
structure functions, and show how the possible 
presence of parton CSV affects those relations. 

Sect. [ITT] reviews both the experimental and the- 
oretical situation regarding CSV in valence quark 
PDFs. In Sect. IIII Al we will summarize recent 
global fits of PDFs that allow for charge symme- 
try violation. Sect. |IIIB| reviews various theoret- 
ical estimates of CSV in valence quark distribu- 
tions. We will argue that one can make reasonably 
model-independent estimates of the magnitude 
and sign of valence parton CSV. In Sect. |IIIB~T] we 
summarize the phenomenon of 'QED splitting,' a 



2 



new source of partonic CSV that results from in- 
clusion of terms where a quark radiates a photon; 
this is the electromagnetic analog of the familiar 
terms where quarks radiate gluons. 



the n th moment of a parton distribution q(x) as 



In Sect. IIIC we summarize the experimental 
limits on valence parton CSV. The most rigor- 
ous upper limits on partonic CSV come from the 
"charge ratio," which compares the F% structure 
functions measured in charge-changing reactions 
induced by neutrinos and antineutrinos, with the 
F2 structure function from charged lepton DIS, in 
principle both measured on isoscalar targets. This 
is reviewed in Sect. Midi In Sect. IIIIC 2l we re- 
view at length the NuTeV measurements of v and 
v reactions on iron targets, and the resulting ex- 
traction of the Weinberg angle. We particularly 
examine potential contributions to this result from 
partonic CSV. 

A new nuclear effect has been proposed that will 
mimic the effects of partonic charge symmetry vio- 
lation. We call this "pseudo CSV," This is defined 
in Sec. Ill D| We show that pseudo CSV effects 
could make significant contributions to analyses of 
the NuTeV experiment. We also discuss how this 
effect could be observed by measuring the nuclear 
dependence of the EMC effect . 

Existing and proposed new experimental facili- 
ties offer several opportunities for dedicated preci- 
sion experiments that could significantly improve 
our chances of observing partonic CSV effects, or 
alternatively of lowering the current upper limits 
on such effects. In Sect. IIII El we summarize four 
such experiments. We show the size of the effects 
that are compatible with the current limits on par- 
tonic CSV, and we discuss those experimental fa- 
cilities that would be best suited to such measure- 
ments. 

In Sect. |IV| we review the situation regard- 
ing sea quark CSV. In contrast to valence quark 
CSV, where there are reliable and rather model- 
independent estimates of the magnitude and sign 
of such effects, it is substantially more difficult ei- 
ther to make theoretical predictions of sea quark 
CSV, or to conceive of experiments to measure 
such effects. In Sect. IIV Al we review theoretical 
and phenomenological estimates of sea quark CSV. 
Perhaps surprisingly, the phenomenological fit by 
Martin, Roberts, Stirling and Thorne (MRST) [23] 
found evidence for a rather large sea quark CSV 
effect. One potentially promising way to test sea 
quark CSV effects involves partonic sum rules. 
The most popular QCD sum rules involve the first 
moment of some combination of structure func- 
tions. For the purposes of this review, we define 



x n 1 q(x)dx 



(5) 



The first moment of valence quark CSV effects is 
necessarily zero, in order to preserve valence quark 
normalization. Thus the only CSV effects to sur- 
vive in this integration are from sea quark CSV. In 



Sect. IV C we review the contributions of partonic 
CSV to the Gottfried, Adler and Gross-Llewellyn 
Smith sum rules. In Sect. HV C 41 we review a new 
sum rule proposed by Ma |32j . Such a sum rule 
would be uniquely sensitive to sea quark CSV ef- 
fects. 

In Sect. [V]we provide a summary and outlook. 



II. RELATIONS BETWEEN HIGH 
ENERGY CROSS SECTIONS AND 
PARTON DISTRIBUTIONS 



General form of high energy cross 
sections 



We can write the cross sections for deep inelastic 
scattering in terms of a set of structure functions, 
which depend on the relativistic kinematics of the 
reaction. Through the quark/parton model, these 
structure functions can in turn be written in terms 
of quark/parton distributions |33j . For simplicity, 
we write the cross sections in leading order (LO) 
in QCD. By now, all phenomenological analyses 
of high energy reactions and structure functions 
work in next to leading order (NLO) or higher 
[3"H [3"5] . At sufficiently high energies, quark mass 
effects are small and can be accounted for with 
quite good precision. Current issues in partonic 
analyses involve data, particularly in neutrino ex- 
periments, in the region where the charm quark 
mass cannot be neglected. One way to deal with 
these issues is to work in a variable flavor number 
scheme (VFNS), where one increases the number 
of active quark flavors at various matching points 
[36] . Recently the CTEQ group has examined the 
effects of quark masses in global fit analyses, par- 
ticularly in extracting the strong coupling constant 
a s from such global fits [37] . The MRST group has 
produced a new set of parton distributions at next 
to next to leading order (NNLO) [34 . In their 
VFNS scheme they introduce discontinuities into 
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their coefficient functions that counter the discon- 
tinuities that arise in their parton distributions at 
the matching points. 

The cross section for scattering of a left (L) or 
right (R) handed charged lepton in neutral current 
(NC) deep inelastic scattering reactions has the 
form 



dx dy Q A 
Mx,y)F](x,Q 2 )] 



\xy 2 Fj{x,Q 2 ) 

Q 2 vi ± at 



(Q 2 + M%) 2sm9 w cos6 v 
x [xy 2 Ff{x,Q 2 ) + f 1 (x,y)Ff(x 1 Q 2 ) 

Q 2 



± h{y)xFf{x,Q 2 )} 



Q 2 + Ml 



2sin0 w cos0 w 
+ h(x, y)F 2 z (x, Q 2 ) ± f 2 (y)xFf(x, Q 2 )] 



(6) 



This process is shown schematically in Fig. [T| In 
Eq. (J©]) a is the electromagnetic coupling, M z is 
the mass of the Z° boson, and 8 W is the Weinberg 
angle. We define the quantities 



fi{x,y) = 1-y 



xyM 2 



Mv) 



y 



(7) 



In Eq. ([7]), M is the nucleon mass. These equations 
are usually evaluated at very high energies where 
xyM 2 « s, so we generally neglect this term; in 
this case for the remainder of this paper we will 
approximate f%(x,y) w fx(y) = 1 - y. 

Either a photon or Z boson can be exchanged in 
this process. The relativistic invariants in Eq. (|6| 
are Q 2 = —q 2 , the square of the four momentum 
transfer for the reaction, x and y. For four momen- 
tum k (p) for the initial state lepton (nucleon), we 
have the relations 



Q 2 



2p.q' 
(k + p) 2 



y 



p • g , 

p ■ k' 



(8) 



Explicit expressions for the structure functions Fi 
in terms of parton distribution functions are given 
in Sec. HTCl below. 
In Eq. (J6), we have 



-1 + 4sin 2 9 V 
4sin^wCos 8 V 




FIG. 1: Schematic picture of deep inelastic scattering 
of charged leptons from a nucleon. Neutral-current 
electroweak interactions involve exchange of a photon 
or Z°. 




FIG. 2: Schematic picture of deep inelastic scattering 
involving the charged-current weak interaction initi- 
ated by charged leptons. An intermediate W is ex- 
changed between the leptons and the nucleon. 



-1 



4 sin w cos 9 W 



(9) 



The most general form of the cross section 
for charged current (CC) interactions initiated by 
charged leptons on nucleons can be written 



d 2 a i :r ] 

dx dy 



us I a 
2 I 9c^ 2 ft 7m2 



2sm z e w (Ml+Q 2 



xy 2 F^ (x,Q 2 ) + h(y)F^ (x,Q 2 ) 



T f 2 (y)xF™ (x,Q 2 ) 



(10) 



This process is shown schematically in Fig. [2] It 
involves a charged virtual W ± of momentum q be- 
ing interchanged between the lepton/neutrino ver- 



v ( v ) 




v ( v ) 



X 



FIG. 3: Schematic picture of deep inelastic scatter- 
ing of neutrinos through neutral-current interactions 
mediated by Z° exchange. 
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tex, and the hadronic vertex. In Eq. (10 1, Mw is 



the mass of the charged weak vector boson. 

Similarly, the cross section for charged current 
interactions initiated by neutrinos or antineutrinos 
on nucleons has the form 



d 2 a^ ] 
dxdy 



2sin 2 9 w {Ml + Q 2 ) 



± Mv)xF™ (x,Q 



(11) 



This process is obtained by interchanging the ini- 
tial and final state leptons in Fig. [2] 

Finally NC reactions initiated by neutrinos or 
antineutrinos have the form 



d 2 cr v{v) 
a u NC 



dxdy "~ V2sin 2 <9 w cos 2 6 W (M 2 + Q 2 ) 
xy 2 F 1 z (x,Q 2 ) + f 1 (y)F 2 z (x 1 Q 2 ) 

± f 2 (y)xF 3 z (x,Q 2 ) . (12) 
This process is shown schematically in Fig. [3] 



Charge Symmetry and Parton 
Distribution Functions 



To obtain the charge symmetry violating parton 
distributions, we introduce the CSV parton distri- 
butions for up and down quarks via 

d n (x) = u p (x)-5u(x); u n (x) = d p (x)-Sd(x) , 

(13) 

and analogous relations for the antiquark distri- 
butions. If the quantities 6u(x) and 5d(x) van- 
ish, then charge symmetry is exact. We assume 
that the strange quark distributions are the same 
in both the proton and neutron, as are the anti- 
strange distributions. There is no theoretical or 
experimental reason to expect strange or charm 
distributions to vary significantly from proton to 
neutron. Except at low x 7 the strange and charm 
distributions are also rather small. 

It is useful to divide light quark parton distribu- 
tions into valence quark and sea quark parts. For 
a given flavor q, the valence quark distributions in 
a nucleon are defined by 



u v (x) 
d v (x) 



u{x) — u{x) 
d{x) — d{x) 



(14) 



We will also use a capital letter to denote the sec- 
ond moment of a parton distribution, i.e. 



U v = 



dxxu v {x) . 



(15) 



The quantity U v in Eq. ( 15 I gives the total momen 



turn carried by up valence quarks in a nucleon. 

For heavy quarks and sea quarks, these often 
appear as linear combinations of the sum and dif- 
ference of quark and antiquark distributions. For 
these we use the notation 



q ± (x) = q{x) ±q(x) 



From Eq. (16 1 it is obvious that u v (x) 



(16) 
-(*). 

However as we will examine in detail features of 
the up and down valence quark PDFs, in this re- 
view we will use the notation of Eq. ( 14 1 for the 
light valence quarks. 

The first moments of the valence quark distri- 
butions obey the quark normalization conditions; 



u^x) dx 



d p (x) dx 



s (x) dx 



d™(x) dx 



u™(x) dx 



c~{x)dx = Q . (17) 



The CSV quantities defined in Eq. ( 13 1 can also be 



decomposed into valence and sea pieces. From the 
definitions of valence quark CSV, and the valence 
quark normalization from Eq. (17 1, it is straight- 



forward to show that the first moment of the va- 
lence quark CSV distributions must vanish, i.e. 



Su v (x)dx— / Sd v (x) dx — . (18) 
(i Jo 



If Eq. ( 18 1 was not true, this would mean that the 



valence quark normalization conditions of Eq. ( 17 1 
could not be satisfied. A consequence of Eq. ( 18 1 
is that 



i f i 

Su(x) dx = Su(x) dx 
i) Jo 



and / Sd(x) dx 
lo 



Sd(x)dx. (19) 



The most precise limits on sea quark CSV are 
derived from QCD sum rules, which involve vari- 
ous moments of the structure functions integrated 
over all x. Some observables will involve the first 
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moment of sea quark parton CSV distributions. 
Thus we can distinguish two different types of par- 
tonic charge symmetry. The first, or "strong form" 
of charge symmetry, is the statement that charge 
symmetry violating parton distributions vanish at 
all x. The "weak form" of charge symmetry cor- 
responds to the assumption that the first moment 
of the CSV sea quark parton distributions is zero, 



5u(x) dx 



8d{x) dx — 



(20) 



even if the parton CSV distributions themselves 
do not necessarily vanish. 

Note that valence quark normalization requires 
that the first moments of heavy quark and an- 
tiquark distributions must be identical. From 
Eq. (17) we see that 



s(x) dx 



s(x) dx 



(21) 



with an analogous relation for charm quarks. This 
simply reflects the statement that the nucleon con- 



tains no net strangeness or charm. From Eq. (21 1 



it is tempting to conclude that the strange quark 
and antiquark distributions should be equal for all 
values of x, e.g., 



s (x) = s(x) — s(x) = , 



(22) 



with an identical relation to Eq. ( 22 1 for the charm 



and anticharm distributions. If all strange quarks 
arise from gluon radiation, Eq. ( |22| would be sat- 
isfied since s and s would always be produced in 
pairs. However, there is no compelling theoretical 
reason why strange quarks cannot arise from other 
sources. For example, in "meson-cloud" models 
UHl US SHI SH H3 S3 Hi], which include virtual 
transitions of a nucleon into a baryon and meson, 
the quark resides in the nucleon and the antiquark 
in the meson. Such models naturally lead to differ- 
ences between quark and antiquark PDFs. We will 
return to this issue in Sect. IIIIC 21 where it will 



be relevant in interpreting the "NuTeV anomaly" 
in the determination of the Weinberg angle. 

Eventually, we would expect partonic charge 
symmetry violation to be calculated directly from 
lattice gauge theory. This would require two addi- 
tional inputs into current lattice calculations. The 
first would be to input different up and down cur- 
rent quark masses. The second will be to include 
electromagnetic interactions into lattice calcula- 
tions. The first of these should be straightfor- 
ward. As far as electromagnetic interactions are 



concerned, there arc presently lattice calculations 
that include electromagnetic effects. For exam- 
ple, Blum et al. estimate light quark masses by 
including electromagnetic interactions and calcu- 
lating pion and kaon mass splittings [45 , using two 
flavors of domain- wall quarks. However, for the 
purpose of testing partonic charge symmetry, one 
must include electromagnetic interactions to suf- 
ficient accuracy to account for all of the major 
effects in the nucleon. A good test would be the 
degree to which the inclusion of electromagnetic ef- 
fects on the lattice can reproduce the experimental 
neutron-proton mass difference of 1.3 MeV. 



Structure Functions in Terms of Parton 
Distribution Functions 



Introducing the CSV parton distributions from 
Eq. (13 1, we can write the leading order expres- 



sions for structure functions without assuming 
charge symmetry. Most tests of charge symme- 
try involve deep inelastic scattering on isoscalar 
targets, which we label as Nq. Such reactions in- 
volve contributions from equal numbers of protons 
and neutrons. So we write the electromagnetic 
and weak structure functions per nucleon on an 
isoscalar target. These expressions are true un- 
der the following conditions. First, we have ne- 
glected contributions from small components of 
the CKM quark mixing matrix |46l [47] . Second, 
we assume that we are working at sufficiently high 
energies and Q 2 such that the quark mass can be 
neglected (this assumption may be inappropriate 
for the charm quark mass, particularly in the case 
of charged-current interactions initiated by neutri- 
nos). These expressions neglect higher- twist con- 
tributions to the structure functions. The PDFs 
depend on the starting scale [i 2 at which they are 
evaluated; in the following equations we do not ex- 
plicitly include the dependence upon the starting 
scale. 

First we provide expressions for the structure 
functions relevant to NC reactions induced by 
charged leptons, given in Eq. 

36F? N °(x, Q 2 ) = 5 [«+(*) + d+(x)} 
+ 2s + (x) + 8c + (x)-ASd + (x)~du + (x). 

(23) 

In the lowest order quark/parton model, the struc- 
ture function is related to the structure func- 



G 



tion F^ p by 



F2 P (x,Q 2 ) 



l+R(x,Q 2 ) 
l + 4M 2 x 2 /Q 



2 2xF?v(x,Q 2 ) 



(24) 

In Eq. (24 1, R = (Jl/ctt is the ratio of the cross 



section for longitudinally to transversely polarized 
photons. An empirical relation fit to the world's 
available data on R has been made by Whitlow et 
al. [10] ■ This fit covers the kinematic region x > 0.1 
and Q 2 < 125 GeV 2 . 

For momentum transfers which are sufficiently 
small (relative to M§) and for parity-conserving 
interactions, we can neglect the contribution from 
Z bosons, in which case the scattering is a function 
only of the two electromagnetic structure func- 
tions, Fj 7 and F 2 , respectively. The cross terms in- 
volving Z bosons are important either at very large 
values of Q 2 , or alternatively for parity- violating 
lepton scattering where the leading terms cancel. 
The structure functions involving photon- Z inter- 
ference have the form 

6F^ N °(x,Q 2 ) = {2g u v -g d v )[u + {x) + d + {x)] 
+ 2g$[2e + (x) - Sd + (x)} - g d v [2s + {x) - 5u + (x)} ; 

2Ff-> Na {x,Q 2 ) = (<# - g d )[u v (x) + d v (x)} 
- gl5d v {x) + g d Ju v (x) . (25) 

The structure functions corresponding to Z° ex- 
change can be written 

AF™>{x, Q 2 ) = (G 2 U + G 2 )[u+{x) + d+(x)] 
+ G 2 d [2s+(x) - Su+(x)} + G 2 u [2c+(x) - Sd+(x)}; 

2F 3 ZN °(x,Q 2 ) = (gl - g d v )[u v (x) + d v (x)} 
+ g^[2c-(x)-Sd v (x)]-g d [2s-(x)-6u v (x)] . 

(26) 



In Eqs. (25 1 and (26 1 we have introduced the re- 
lations 



2 " 3 Sm ° w 

2 ■ 2n 1 



u\2 



Gl = (g 
G 2 d = (gt) 2 + 



1 



(27) 



For charged-current interactions initiated by ei- 
ther charged leptons or by neutrinos, for suffi- 
ciently low values of Q 2 one must take account 
of the masses of heavy quarks. We do not in- 
clude bottom and top quark effects in this re- 
view; however one must account in some way for 
the nonzero charm quark mass. In this case the 



charged-current structure functions will depend on 
the CKM matrix elements [1HJ 27] ■ Expressions 
for the charged-current structure functions that 
take into account effects of heavy quark masses 
and CKM matrix elements can be found in the lit- 
erature [33]. However, for sufficiently lar ge va lues 
of Q 2 , the structure functions in Eqs. ( 10 ) and (111 



will to a good approximation simplify to the form 

2F^ +No (x, Q 2 ) -> u+{x) + d+(x) + 2s(x) 
+ 2c{x) — 5u{x) ~ 8d(x) ; 

2F^ Na {x, Q 2 ) -> u+(x) + d+(x) + 2s(x) 
+ 2c(x) — 5d(x) — 5u{x) ; 

FY +No {x, Q 2 ) -> u v (x) + d v (x) + 2s{x) 
— 2c{x) — 5u{x) + 5d(x) ; 



Ff No (x,Q 2 ) -» u v (x) + d v (x) - 2s(x) 
+ 2c(x) - Sd(x) + Su(x) . 



(28) 



In Eq. (23 1 and subsequent equations, we have 



suppressed the nucleon index on the parton distri- 
butions. Since we have explicitly introduced the 
parton CSV amplitudes, the remaining PDFs are 
now understood to be those for the proton. The 
relation between the F\ and F% structure functions 
for neutrinos is given by an equation analogous to 
Eq. (24 1, where R is now the longitudinal to trans- 
verse ratio that holds for CC and NC neutrino 
reactions. For the CC reactions initiated by neu- 
trinos, the experimental values for R v are summa- 
rized in Conrad, Shaevitz and Bolton [48]. For NC 
reactions, the value of R is essentially unknown. 
This provides some uncertainty in extracting par- 
ton distribution functions from neutrino NC reac- 
tions [4*5] . 

We have written the nuclear structure functions 
in terms of the parton distributions for free nucle- 
ons. However, as is well known, parton distribu- 
tions are modified in nuclei. At small x there are 
shadowing corrections, at intermediate x there are 
'EMC effects,' [5H5TJ [52] , and at large x Fermi- 
motion effects dominate. Nuclear modifications of 
PDFs have been reviewed recently by Kumano and 
collaborators [53] [54] [55] , and also by Kulagin and 
Petti [55] [57] [55] . Consequently, in any precision 
experiments these effects must be accounted for 
if we compare to parton distributions taken from 
free protons (there should even be small modifica- 
tions arising in the deuteron [5"3 [50] 1E5 ) • We will 
discuss these effects later as they arise. 
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D. Relations Between Structure Functions 

Using the relation between leading-order high- 
energy structure functions given in Eqs. (23 1 and 



( 28 1 , we obtain the following relation between the 



structure functions, including charge symmetry vi- 
olating effects. 



18 

xF. 



-F 



WN 0l 



{x)-Fl No (x) = ^\xF™ +N °(x) 



W~N 



(X) 



X 

12 



2(s+(x)-c+(x)) 



+ 6d + (x)-5u+(x) 



(29) 



Eq. ([29} , sometimes called the "5 /18 th rule," 
relates the F2 structure function from charged- 
current neutrino reactions to the F 2 structure 
function from interactions of charged leptons, with 
both quantities measured on isoscalar targets. In 
Eq. ( 29 ) we have for simplicity neglected the lon- 



gitudinal to transverse correction factors R given 
in Eq. ( 24 1 . However these correction factors were 



included when the F2 structure functions were ex- 
tracted from experimental cross sections. 
1 1 — WN 

In Eq. ( 29 1, F 2 is the average of the CC cross 



sections induced by v and v, 



F 



WN 



F 2 w+No {x) + F 2 w ~ No (x) 



(30) 



The right-hand side of Eq. ( 29 1 includes contribu- 



tions from strange and charmed quarks, and is cor- 
rect to lowest order in CSV terms. Although the 
light quark contributions cancel in this expression, 
there is a residual contribution from strange and 
charm quarks. At small x one has contributions 
from both the CSV and heavy quark PDFs. How- 
ever, at larger values of x the strange and charm 
contributions should be extremely small, and in 
this region the only significant contribution to the 
right-hand side of Eq. ([29l should come from va- 



lence quark CSV terms. Furthermore, if the heavy 
quark contributions are known, Eq. (29 1 may be 



used to investigate the charge-symmetry violating 
quark distributions for the light quarks. 

The strange quark PDFs have been determined 
from the production of opposite-charge muon pairs 
in neutrino-induced reactions [62J, [63l [64j [65j |66l 
|571|55]. Thus comparison of these two F 2 structure 
functions, combined with our knowledge of strange 
and charm PDFs, has the potential to measure (or 
to place strong upper limits on) parton CSV prob- 
abilities. The current experimental and theoreti- 
cal situation will be reviewed in Sec. IIII C II From 



Eq. ( 28 1 , we note that in principle we could obtain 
the same information as in Eq. ( 29 1 , by measur- 



ing the difference between the xF$ structure func- 
tions from charge-changing v and V interactions 
on isoscalar targets. 

We can obtain another relation between struc- 
ture functions by measuring the F 2 structure func- 
tion from v and V CC reactions on isoscalar tar- 
gets. Using Eq. (28 1 we obtain 



F. 



W+N 



2 

x[2(s" 



F 2 "°(z) 



(x) 

(x) — c~(x)) + Sd v (x) — Su v (x)]. 



(31) 



The right-hand side of Eq. (31) contains "va- 
lence" contributions (the difference between quark 
and antiquark probabilities) for strange and charm 
quarks, as well as contributions from valence quark 
CSV terms. In Sec. Ill E 4| we discuss the ex- 



perimental possibilities for measuring this quan- 
tity, and we show theoretical predictions and phe- 
nomenological limits on this quantity. 

Another relation for structure functions can be 
obtained by comparing the F 2 structure functions 
obtained from charge-changing interactions of an- 
tineutrinos and neutrinos on proton targets. If one 
takes the difference between these structure func- 
tions and divides by 2x one obtains 



F, 



W~p 



(x)-F 2 w P(x) 



2x 

u v (x) — d v (x) — s~(x) + c~{x). (32) 

Since the first moment of the strange and charmed 
"valence" contributions must vanish, the differ- 
ence between these two F 2 structure functions, di- 
vided by 2x and integrated over all x, should be 
equal to the difference between the up and down 
valence quark occupation numbers in the proton, 
or one. This is the Adler sum rule [69], which will 
be reviewed in Sec. IIV C 21 

If the -F3 structure functions for neutrino and 
antineutrino charge-changing reactions are mea- 
sured on isoscalar targets, then from Eq. (28 1 the 
sum of these structure functions gives 



xF^ +N °(x) + xF^ N °(x) 



2x 

+ s~(x) + c~(x) - 



u v (x) + d v (x) 



Sd v (x) + 5u v (x) 



(33) 



The sum of these structure functions includes only 
valence quark probabilities plus valence CSV con- 



tributions. Consequently integrating Eq. (33 1 over 
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all x, and applying valence quark normalization 
from Eqs. (17 1 and (18 1 gives (modulo QCD cor- 
rections) just the sum of valence up and down 
probabilities in the nucleon, or three. This is the 
Gross-Llewellyn Smith sum rule |70j . reviewed in 
Sec. HSTH 

One final relation can be obtained by comparing 
the F2 structure function from charged lepton DIS 
on protons with that for neutrons. One obtains 



+ 



F 2 7P (x) - F2 n {x) _ u+{x)-d+{x) 

x 3 
45d+(x) + 5u+(x) 



(34) 



If the quantity in Eq. (34 1 is integrated over all x 



then one obtains the Gottfried sum rule [TTj. The 
experimental and theoretical implications of this 
sum rule are discussed in Sec. II V C~T1 



III. CHARGE SYMMETRY IN VALENCE 
QUARK DISTRIBUTIONS 



In this section, we will review both theory 
and experiment regarding charge symmetry in va- 
lence quark distributions. First, we will review 
phcnomcnological estimates of CSV for valence 
quarks. Next, we will review theoretical estimates 
of valence CSV, and then we will review the exper- 
imental upper limits. Finally we will suggest new 
experiments that could provide strong constraints 
on valence quark CSV. 



A. Phenomenological Estimates of Valence 
Quark CSV 

Recently, Martin, Roberts, Stirling and Thorne 
(MRST) [23 j have evaluated uncertainties in par- 
ton distributions arising from a number of factors, 
including isospin violation. They chose a specific 
model for valence quark charge symmetry violat- 
ing PDFs, adopting a function of the form 



Su v (x) 



-Sd v (x) 



(1 



\4 - 
X) X 



-- Kf(x) ; 

" 5 (a - .0909) 



(35) 



The quantity f(x) in Eq. (35) was chosen so that 
its x dependence had roughly the same form as 
the MRST valence quark parton distribution func- 
tions (at the starting scale for QCD evolution) in 
both the limits x — > and x — > 1. The first mo- 
ment of f(x) was fixed to be zero, in agreement 



with the valence quark normalization constraint 
of Eqs. (17 1 and (18 1. The valence quark normal- 



ization condition requires that the CSV function 
f(x) have at least one node. 

Inclusion of valence quark CSV can in principle 
change the momentum carried by valence quarks 
in the neutron from those in the proton, since 
the total momentum carried by valence quarks is 
given by the second moment of the distribution. 
The total momentum carried by valence (up plus 
down) quarks in the neutron is determined exper- 
imentally to within about 2%, so MRST chose a 
functional form that insured that equal momen- 
tum was carried by valence quarks in the proton 
and neutron. For this reason, they insisted that 
the valence CSV terms Sd v and Su v be equal and 
opposite. With this constraint it is straightfor- 
ward to show that 



u v v (x) + d%{x) = v%(x) + d%{x) 



(36) 



Multiplying both sides of Eq. ( 36 1 by x shows that 



the momentum carried by valence quarks is iden- 
tical for proton and neutron. The overall coeffi- 
cient k in Eq. ( 35 1 was varied in a global fit to 



a wide range of high energy data. For simplicity, 
MRST neglected the Q 2 dependence of the CSV 
term in their global fit. Later, when we use these 
charge symmetry-violating PDFs to estimate po- 
tential effects of partonic CSV, this will introduce 
some uncertainty. To lowest order in QCD, typ- 
ical CSV effects will have a form proportional to 
something like 



Su(x) — 5d(x) 
u(x) + d(x) 



(37) 



If we use the MRST CSV PDFs obtained using 
Eq. (35 1, for a given Q 2 we will be using par- 



ton distributions where the denominator has been 
evolved in Q 2 but the numerator is not evolved. 

Including valence quark CSV in their global fit 
to high energy data, MRST obtained a very shal- 
low minimum in \ 2 with a best-fit value k = —0.2, 
and a 90% confidence level for the range —0.8 < 
k < +0.65. The x 2 for their fit vs. the parameter 
k is shown in Fig. |4j Since MRST chose a very 
specific functional form for valence quark CSV, 
their results could have a substantial model de- 
pendence. The MRST global fit guarantees that 
CSV distributions with this shape, and with val- 
ues of k within the 90% confidence range, will give 
reasonable agreement with all of the high energy 
data used to extract quark and gluon PDFs. 

Since the MRST functional form for valence 
CSV PDFs requires that 5d v be equal in magni- 



9 




FIG. 4: [color online] The \ 2 obtained by MRST, Ref. 
|23| for a global fit to high energy data of parton dis- 
tribution functions including valenc e q uark CSV with 
the functional form defined in Eq. ( 35 1 
vs. the free parameter re. 



X 2 is plotted 



tude to 6u v , this implies that at large x the frac- 
tional charge symmetry violation is substantially 
larger for the "minority valence quark" distribu- 
tion d v than for u v , since d v << u v in this region. 
Similar results have been obtained for valence CSV 
distributions within a number of theoretical mod- 
els, as will be discussed in the following section. 
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FIG. 5: The phenomenological valence quark CSV 
function from Ref. [22] , corresponding to best fit value 
re = —0.2 defined in Eq. (35 1. Solid curve: xSd v ; 
dashed curve: xSu v . 

In Fig. [5] we show the phenomenological valence 
quark CSV distributions obtained by MRST us- 
ing the function f(x) from Eq. (35 1 with the pa- 
rameter k = —0.2, which represents the best fit 
to the high energy data. The solid curve cor- 
responds to xSdyQr) and the dashed curve cor- 
responds to xSu v (x). These valence CSV PDFs 



reach a maximum value of approximately 0.006 at 
a value x ~ 0.3, and (by inspection of Eq. (35 1) 
they have a zero crossing at x = 0.0909. 



B. Theoretical Estimates of Valence Quark 
CSV 

The phenomenological MRST results of the pre- 
ceding section can be compared with theoretical 
estimates of valence quark CSV. In a valence quark 
approximation, the nucleon can be considered as 
consisting of three valence quarks, with proton and 
neutron described as 



\p) 
In) 



[uud] , 
[udd] . 



(38) 



In quark models evaluated on the light cone, the 
valence quark distribution can be expressed as |721 

Eg 



X 



x 5(M(l-x) -Px) ■ 



(39) 



In Eq. (39 1 a 3 = 77; this relation denotes the 
process where a valence quark is removed from a 
nucleon |iV), and the result is summed over all 
final states \X). The quantity p x is the energy 
of the state following removal of a valence quark 
with momentum k. The quantity /1 2 represents the 
starting value for the Q 2 evolution of the parton 
distribution. Eq. (39) treats only the quark longi- 



tudinal momentum and neglects transverse quark 
momentum. 

There are several potential sources of charge 
symmetry violation in Eq. (39). First, there is 



possible charge symmetry violation in the quark 
wavefunctions. Second, there are mass differences 
in the spectator multiquark system. And finally 
there are additional electromagnetic effects which 
break charge symmetry. Now, model quark wave- 
functions are found to be almost invariant under 
the small mass changes typical of CSV [73] , so we 
do not discuss these effects further. Electromag- 
netic effects are of order a where a is the electro- 
magnetic coupling constant; hence such effects are 
expected to be at the 1% level. At the large val- 
ues of Q 2 characteristic of high-energy reactions, 
typically one to many GeV 2 , such effects should 
be small. 

For simplicity, one can consider charge symme- 
try violation to arise from two effects. The first is 
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the n — p mass difference SM = M n — M p = 1.3 
MeV. This mass difference is the result of a num- 
ber of electromagnetic effects, in addition to the 
u — d quark mass difference. A second effect giv- 
ing rise to charge symmetry violation is differences 
in diquark masses arising from the current quark 
mass difference between up and down quarks. We 
define the quantity 



Sin = m dd - m uu 



(40) 



One has a robust estimate for this mass difference, 
Sfh ~ 4 MeV [75] • With these approximations 
one can write charge symmetry violating valence 
parton distributions in terms of 



Sq v 



dq v 
dSm 



5m + 



dg v 
dM 



SM 



(41) 



From Eq. (41) the valence charge symmetry vi- 
olating parton distributions are obtained by tak- 
ing variations with respect to diquark and nucleon 
masses on valence parton distributions from quark 
models. 

If we take the simple valence quark picture of the 
nucleon as given by Eq. ( 38 1 , then we can consider 



diquark mass differences following the removal of 
one quark from the nucleon. If we remove a "ma- 
jority" valence quark (a u quark in the proton or 
d quark in the neutron), then for both proton and 
neutron one is left with a ud diquark. Thus for 
the majority quark distribution, there is no quark 
mass asymmetry for the residual diquark. For re- 
moval of a "minority" quark, (a d quark in the 
proton or u quark in the neutron), the remainder 
is a uu diquark in the proton and a dd diquark in 
the neutron. Thus the diquark mass asymmetry is 
just given by the quantity Sfh in Eq. ( 40 1 . 

This technique was used by Sathcr |76J, who in- 
vestigated these effects in a static quark model. 
Sathcr obtained an analytic approximation relat- 
ing valence quark CSV to derivatives of the valence 
PDFs 

. . . SM ( d ... d , . 
ou v {x) = — I -— [xu v (x)\ + ^M^) 



(42) 

Note that Sather's equations agree with our earlier 
arguments. The "majority" quark CSV distribu- 
tions Su v (x) — v%(x) — rf"(ar) are functions only of 
SM and do not depend on Sfh, while the "minor- 
ity" valence quark CSV distributions 5d v (x) de- 
pend on both SM and Sfh. 



In Fig. [6] we plot the CSV valence parton distri- 
butions of Sather (75] . The valence parton distri- 
butions used by Sather were from preliminary fits 
to Tevatron structure functions obtained by the 
CCFR group. The solid curve is xSd v (x) vs. x, 
while the dashed curve plots xSu v (x). The PDFs 
have been evolved to Q 2 = 12.6 GeV 2 . The quali- 
tative features of these charge symmetry violating 
valence PDFs are similar for all models that we 
will review. Since the first moment of the valence 
CSV parton distributions has to vanish (in order 
to maintain valence quark normalization), the va- 
lence CSV PDFs must change sign at least once. 
For Sather's model this occurs at x ~ 0.05. In 
general, if one calculates the CSV valence PDFs 
by inserting phenomenonogical parton distribu- 
tion functions into Eq. (42 I or other analytic for- 



mulae, the resulting CSV parton distributions will 
not obey the valence quark normalization condi- 
tion. However, Sather obtained his PDFs from 
moments of the quark distributions; for his va- 
lence CSV parton distributions Sather set the first 
moment to zero, thus guaranteeing valence quark 
normalization. 
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FIG. 6: Theoretical valence CSV PDFs from Sather, 
Ref. 76 and Eq. (|42]). Solid curve: x5d v ; dashed 
curve: x5u v . The valence PDFs were preliminary 
fits to CCFR Tevatron structure functions evolved to 
Q 2 = 12.6 GeV 2 . 

For larger values of a;, Sd v (x) is positive while 
Su v (x) is negative. The distributions peak at 
x ~ 0.3. In Sather's model Sd v (x) is roughly 50% 
larger in magnitude than Su v (x). By observing 
the quark model wavefunctions we can understand 
these qualitative features of the valence quark CSV 
distributions. The "minority valence quark" CSV 
term is defined by Sd v = d^ — u^. Removing a mi- 
nority valence quark from a nucleon with 3 valence 
quarks leaves a diquark system that is uu for the 
proton and dd for the neutron. Simple theoretical 
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arguments [71] suggest that the down quark dis- 
tribution in the proton will be shifted to higher x 
and the up quark distribution in the neutron will 
be shifted to lower x. This predicts that, at large 
x, 5d v should be positive. 

Conversely for the "majority valence quark" 
CSV term 5u v = v% — d™, removing a majority 
quark leaves intermediate states with the same 
quark configuration (ud) for both neutron and pro- 
ton. From Eq. (42 I the majority valence CSV term 
should depend only on the n — p mass difference, 
and one expects that 5u v should be negative at 
large x. These qualitative predictions agree with 
the quark model CSV valence distributions shown 
in Fig. [6] and they also are in agreement with 
the phenomcnological best-fit CSV PDFs shown 
in Fig. [5] Although the down valence distribu- 
tion in the proton is less than half the up valence 
distribution, these qualitative arguments suggest 
that 



Sd v (x) > \Su v (x)\ 



(43) 



at large x. This is observed in the theoretical 
model calculations by Sather, however Eq. (43 1 is 



not satisfied by the phenomenological MRST pa- 
rameterization of Eq. ( 35 1 , which requires by def- 
inition that Sd v (x) and Su v (x) should be equal in 
magnitude and opposite in sign. 

There is another way to understand the quali- 
tative features of these valence CSV distributions. 
It was pointed out by Londergan and Thomas |77j 
that from Sather 's expression Eq. (42 1 one can ob- 



tain an analytic expression for the second moment 
of the CSV parton distributions, 



6U V 
6D V 



SM 
SM 



(U v - 2) , 
Sfh 



D 

M M 
SM 



M 



(D v + 3) 



(44) 



In Eq. (44), U v and D v are respectively the to- 
tal momentum carried by up and down valence 
quarks. The final line of Eq. ([44]) follows from the 
fact that Sfh ~ 3SM. Since the valence CSV dis- 
tributions are required by valence quark normal- 
ization to have zero first moment (see Eq. ( 18)), 



the valence CSV distributions must change sign at 
least once. Eq. (44 1 predicts that at large x 5u v 



will be negative and Sd v will be positive. The sign 
of the valence CSV distributions is the same for 
all parton distributions derived from quark mod- 
els. Furthermore, from the relative magnitude of 



SU V and SD V , we expect the maximum of Sd v to 
be larger than that for Su v . 

Benesh and Londergan [78] also considered par- 
ton charge symmetry violation from quark mod- 
els using Eq. (39 1. They related the change in 



PDFs due to charge symmetry violation in minor- 
ity valence quark distributions (the term propor- 
tional to the diquark mass difference Sfh), to the 
color hyperfine splitting between N and A states 
in quark models of baryons that initially assume 
SU{A) spin-flavor symmetry (these arise from di- 
quark spin splittings). Using the work of Close 
and Thomas [75], they obtained 



Sd v (x) 
Su v (x) 



5m 
S h f 



u v (x) — 2d v (x) 



SM d 
M dx 



G 



U V (X) 



SM d 
M dx 



dy{x) , 



(45) 



In Eq. (|45|), S hf = 50 MeV is the S = 1 color hy- 
perfine splitting in the SU(4) limit. These equa- 
tions also differ from Sather's result of Eq. (42 1 in 



that Benesh and Londergan considered variations 
of the nucleon mass M while keeping the quantity 
Mx constant, following arguments by Benesh and 
Goldman [50]. 

In Fig. [JJ we plot the theoretical valence CSV 
parton distributions calculated by Benesh and 
Londergan [78], using Eq. pS). The PDFs used 



were the phenomenological CTEQ LQ (low Q) 
parton distributions from the CTEQ group [HI], 
evaluated at the low momentum starting scale 
Q 2 = 0.49 GeV 2 . 

As we have mentioned (see Eq. ( [181 ), va lence 
parton CSV distributions should respect valence 
quark normalization, and hence (Sd v ) = (5u v ) = 
0. By inspection of Fig. [7] the valence CSV PDFs 
of Benesh and Londergan do not satisfy the quark 
normalization condition. Although the sign of the 
CSV PDFs obtained by Benesh and Londergan 
agrees with that of Sather, and also with the pre- 
dictions of Eq. ( 44 1 , the magnitudes are somewhat 
different. Here the magnitude of Su v is larger than 
Sd v , a result opposite from Sather and in disagree- 
ment with Eq. (44 1. The difference between these 



theoretical PDFs is likely related to the fact that 
Benesh and Londergan used an additional approx- 
imation to relate quark CSV terms to mass split- 
tings in SU(4) symmetric quark models. Note that 
the Benesh-Londergan PDFs are evaluated at a 
considerably lower value of Q 2 than for Sather. 
One could evolve these CSV PDFs to higher Q 2 by 
inserting these parton distributions into the QCD 
evolution equations. 
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FIG. 7: Theoretical valence CSV PDFs from Benesh 
and Londergan, Ref. [78] and Eq. (45 1. Thick solid 
line: xSd v ; thin solid line: x8u v . Dash-dot line: xSd v 
when the n-p mass difference is set to zero. The PDFs 
were the CTEQ LQ (low Q) parton distributions from 
Ref. |81| . evaluated at the low momentum starting 
scale Q 2 = 0.49 GeV 2 . 



Rodionov, Thomas and Londergan [71] also cal- 
culated charge symmetry violating parton distri- 
butions using Eq. (39 1. They also included the 



quark relativistic energy (if one ignores this the 
resulting parton distributions do not have the cor- 
rect support; the PDFs are then defined on the 
range < x < oo, rather than from < x < 1). 



Rodionov et al. also evaluated Eq. (39 1 including 
the effects of quark transverse momentum. In this 
case one can no longer obtain analytic expressions 
for the CSV valence parton distributions. 
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FIG. 8: Theoretical CSV PDFs by Rodionov et al, 
Ref. [74 . Solid line: x8u v ; dash-dot line: xSd v . The 
PDFs have been evolved to Q 2 = 10 GeV 2 . 

Fig. |8]plots the theoretical quark model calcula- 
tions of valence CSV by Rodionov et «Z. [75]- The 



dash-dot curve in Fig. [8] represents the quantity 
xSd v (x) while the solid curve is xSu v (x). The 
curves were initially calculated at the low Q 2 ap- 
propriate for quark models, and evaluated at Q 2 = 
10 GeV 2 through DGLAP evolution [55] [57J [2B] . 
The valence CSV parton distributions obtained by 
Rodionov et al. |74] are quite similar to those of 
Sather [75], as seen by comparison of Figs. [6] and 
§ The sign and magnitude of both Sd v (x) and 
5u v (x) are very similar, and the second moments 
of both distributions agree to within about 20%. 
The zero crossing in Sather's model appears at a 
smaller value of x than that for Rodionov. 

We can also compare the theoretical valence 
CSV distributions with the phenomenological va- 
lence CSV distributions obtained by MRST from 
their global fit to high energy data. These are 
plotted in Fig. [5] for the best fit value k — —0.2 
in Eq. (351 . The solid (dashed) curve in Fig. [5] 
represents x5d v (x) (xSu v (x)). The sign and rel- 
ative magnitude of both 6d v and 5u v , and the 
point where they cross zero, are remarkably similar 
in both the MRST phenomenological CSV PDFs, 
and the results obtained by Rodionov et al. and 
shown in Fig. [8] The second moments of both of 
the Rodionov quark model CSV PDFs are equal 
to the moments of the corresponding MRST values 
to better than 10%. 

The valence CSV parton distributions Su v ob- 
tained by Benesh and Londergan [78j and shown 
in Fig. [7J are quite similar to those of Sather and 
Rodionov, while the CSV valence distribution Sd v 
is roughly a factor of two smaller than the others. 
Benesh and Goldman [HO] calculated parton CSV 
distributions from a quark potential model, and 
their CSV PDFs have the same sign and a similar 
shape to those derived by Sather and Rodionov, 
but the Benesh- Goldman CSV PDFs are roughly 
a factor two smaller in magnitude. 

The qualitative agreement between the phe- 
nomenological valence quark PDFs obtained by 
MRST, using the best value k — —0.2 from their 
global fit and shown in Fig. [3J and the theoretical 
CSV PDFs obtained by Rodionov et al. and shown 
in Fig. [8] is rather remarkable, especially consider- 
ing that the theoretical results were obtained some 
ten years earlier and used relatively simple bag 
model quark wavefunctions. The excellent agree- 
ment with the phenomenological results provides 
some theoretical support for the functional form 
chosen by MRST. However, within the 90% confi- 
dence region for the global fit, the valence quark 
CSV PDFs could be either four times as large as 
those predicted by Sather and Rodionov, or they 
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could be three times as large with the opposite 
sign. 

One feature of the theoretical CSV distributions 
is the prediction that for moderately large values 
of x (i.e. for x above the zero crossings) 

\5u v (x) + Sdy(x)\ « \5u v (x) - 8dy(x)\ . (46) 

Consequently, valence quark CSV observables that 
depend upon the difference between the minor- 
ity and majority CSV terms should be substan- 
tially larger than those that depend on the sum 
of these terms. Eq. (46 1 is satisfied trivially for 



the MRST phcnomcnological valence CSV distri- 
butions of Eq. ( 35 1 , since by definition their sum 
is zero. 

Cao and Signal [82 calculated partonic charge 
symmetry violation assuming that partonic CSV 
arises through mesonic fluctuations of the nucleoli. 
They used a meson-cloud model to estimate par- 
tonic CSV [HI |53]- In the meson-cloud model, 
mass splittings in the baryon and meson multiplets 
lead to charge symmetry violating parton distribu- 
tions. The resulting CSV distributions obtained 
by Cao and Signal are substantially smaller than 
those obtained by Sather [7S] or Rodionov [73], 
and peak at substantially smaller values x ~ 0.1. 
This could be expected from the splitting functions 
for baryons in meson-cloud models. 

Cao and Signal break up the parton distribu- 
tions into three parts: a bare part, a perturba- 
tive part and a non-perturbative part. The first 
two of these are assumed to be charge symmetric. 
The perturbative part is assumed to arise from 
gluon splitting. In the next section we will show 
that there is an additional perturbative part aris- 
ing from photon splitting. This additional part 
will contribute to parton charge symmetry viola- 
tion since the photons couple differently to up and 
down quarks, by virtue of their different charges. 



"QED Splitting:" Another Source of Parton CSV 



Recently, another source of parton charge sym- 
metry violation has been included in calculations 
of PDFs by both MRST [24] and Gliick, Jimenez- 
Delgado and Reya [25] . The most important terms 
in the usual QCD evolution involve gluon radi- 
ation, where a quark radiates a gluon leaving a 
quark with a lower x value. These authors sug- 
gested that one assume charge symmetry at some 
initial extremely low-mass scale, and include in the 



QCD evolution equations the effect of photon radi- 
ation. This involves the explicit coupling of quarks 
to photons, the analog of quark coupling to gluons. 

Fig. [9] presents a schematic picture of the cou- 
pling of quarks to photons. The electromagnetic 
couplings are obtained by replacing gluon lines 
with photons, except for the gluon self-coupling 
which is not present for photons, and replacing the 
gluon splitting functions with the appropriate cou- 
pling for photons. This QED coupling changes the 
parton distribution functions in two distinct ways. 
First, it introduces an additional source of charge 
symmetry violation, since the photon couples dif- 
ferently to up and down partons because of their 
different electromagnetic charges. Second, radia- 
tion of the photon produces a "photon parton dis- 
tribution." This photon PDF must be accounted 
for in the evolution equations. The photon PDF 
also makes a contribution to the total momentum 
carried by the nuclcon. 
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FIG. 9: Schematic picture of quarks coupling to pho- 
tons. Replacing a gluon line by a photon everywhere 
(except for the gluon self-coupling) produces the elec- 
tromagnetic coupling of photons to partons. This gives 
the origin of QED splitting that produces additional 
CSV effects in parton distribution functions. 

When one includes QED contributions in this 
way, to lowest order in both the strong coupling 
a s and the electromagnetic coupling a, the so- 
called DGLAP evolution equations due to Dok- 
shitzer [53], Gribov and Lipatov [27] and Altarelli 
and Parisi [25] are modified. To lowest order in the 
QED coupling a the evolution equations obtained 
by MRST have the form 
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The convolution integral in Eq. [47] is denned by 
rl dy 



P®f 



(48) 



In Eq. ( 47 1 , the right hand side of the schematic 



evolution equations represents a convolution of the 
splitting functions with the quark and gluon dis- 
tributions (which have an explicit dependence on 
the factorization scale parameter ji 2 ). Inclusion of 
the electromagnetic contribution to the evolution 
equations introduces a "photon parton distribu- 
tion" 7(x,^ 2 ) which is coupled to the quark and 
gluon distributions. The new splitting functions 



that occur in Eq. ( 47 1 are related to the standard 



QCD splitting functions by 
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(49) 



In Eq. (49 1, N c is the number of colors. Conser- 



vation of momentum is assured by the relation 



dxx 



1 



It is necessary to simplify Eq. (47) 



(50) 
First, since 

the electromagnetic interaction is not asymptot- 
ically free, it is difficult to determine a model- 
independent method for setting the starting val- 
ues for the various PDFs that are coupled by these 
QED effects. In particular, it is not clear where the 
QED effects should be assumed to vanish. Second, 
inclusion of the QED couplings could in principle 
more than double the number of parton distribu- 
tion functions (one must now differentiate between 
proton and neutron PDFs, in addition to the new 
photon PDFs). 

Two groups have adopted somewhat different 
strategies, with similar overall results. Gliick et 
al. [25] adopt the standard convention for DIS re- 
actions of setting the scale fi 2 = Q 2 . The most 
important contribution from the photon coupling 
occurs in the valence quark PDFs. To lowest order 
in the electromagnetic coupling a, the convolution 
equations for the CSV valence quark distributions 
arising from QED coupling have the form 
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(51) 



Similar relations hold for the antiquark distribu- 
tions. Gliick et al. assume that the average cur- 
rent quark mass WL q = 10 MeV is the kinematical 
lower bound for a quark to emit a photon. This is 
analogous to taking the electron mass as the lower 
limit for radiation of photons in the earliest cal- 
culations of the Lamb shift (before the advent of 
renormalization group arguments) |84j . Eq. (51 1 is 



then integrated from m 2 to Q 2 . QED evolution ef- 
fects are evaluated while keeping the QCD effects 
fixed. The quark distributions appearing on the 
right hand side of Eq. (51 1 are the GRV leading- 
order parton distributions 85j. In the resulting 
integrals, in the region m 2 < q 2 < fi 2 LO = 0.26 
GeV 2 corresponding to momentum transfers below 
the input scale for GRV, the PDFs are "frozen," 
i.e. in this region they are assumed to be equal to 
their value at the input scale ^\q- 
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FIG. 10: The isospin-violating majority xSu v (solid 
curve) and minority xSd v (dashed curve) valence par- 
ton distributions obtained by Gliick et al. [25] at Q 2 = 
10 GeV 2 , assuming QED evolution from a scale set by 
the current quark mass. These are compared with ma- 
jority (solid points) and minority (open circles) CSV 
distributions obtained from theoretical quark model 
calculations [74] . 

The resulting valence isospin asymmetries xSu v 
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and xSd v are plotted in Fig. 10 at Q 2 = 10 
GeV 2 . For comparison, they are plotted along 
with the valence quark CSV asymmetries obtained 
from quark model calculations by Rodionov et 
al. |711 US]. In the Rodionov calculations CSV 
distributions arose from diquark mass differences 
Sfh = nidd — m uu ss 4 MeV [TH [76] , and from the 
target nucleon mass difference 5M = M n — M p . 
While the quantity Su v is quite similar in both 
sign and magnitude for both the bag model and 
the QED calculations, the QED results for Sd v are 
roughly half as large as the bag model results. This 
can be understood from the evolution equations of 
Eq. (51). The coefficients of the QED evolution 



are equal and opposite for up and down valence 
quarks, but since u v is roughly twice d v , one ex- 
pects the CSV effects for up quarks to be approxi- 
mately twice the magnitude and the opposite sign 
as those for down quarks. As a result, the CSV ef- 
fects obtained from "QED splitting" will not obey 
the relation of Eq. ( 36 I assumed by MRST in their 



phcnomenological fit. As noted previously, the bag 
model results for valence quark CSV are extremely 
close to those obtained by MRST using the phe- 
nomenological form of Eq. (35), for the best-fit 
value k — —0.2. 

The MRST group [H] solved the evolution equa- 
tions of Eq. ( 47 ) with assumptions about the par- 



ton distributions at the starting scale Qq = 1 
GeV 2 . At the starting scale, the sea quark and 
gluon distributions are assumed to be charge sym- 
metric. The photon PDFs at the starting scale 
were taken as those due to one-photon radia- 
tion from valence quarks in leading-logarithm ap- 
proximation, evolved from current quark masses 
vn u — 6 MeV and m d = 10 MeV to Qq. This 
produces different photon PDFs for the neutron 
and proton at the starting scale. Enforcing over- 
all valence parton momentum conservation from 



Eq. ( 50 ) requires valence quark charge asymme- 



try at the starting scale. MRST assumed a simple 
phcnomenological form chosen to obey the valence 
quark normalization condition, which produced 
charge symmetry violating distributions that re- 
semble the valence PDFs at large and small x, 
and with an overall magnitude chosen to enforce 
quark momentum conservation. MRST then de- 
termine the proton's quark and gluon distributions 
at the starting scale Qq by a global fit to an ar- 
ray of high energy data. The only difference from 
other MRST global fits is the use of the modi- 
fied DGLAP evolution equations of Eq. ((47]). The 



The contribution to CSV arising from QED 
splitting would occur even if the up and down 
quark masses, and the neutron-proton masses, 
were initially identical. This is different from the 
CSV terms which were calculated from quark mod- 
els and described in the preceding section. From 
Eq. (41 1, it is clear that those CSV terms were 



proportional to the up-down quark mass difference 
and the n — p mass difference. 

Because the two types of parton charge symme- 
try violation tend to arise from different sources, 
and both CSV effects are quite small, we have 
evaluated them independently and we add them 
together. Note, however, that the CSV contri- 
butions from QED splitting cannot be treated as 
being completely independent of the quark-model 
CSV terms. This is because the quark-model cal- 
culations used estimates of electromagnetic effects 
in calculating mad — fn%u a s described following 



Eq. (421 



The quark PDFs calculated using the QED 
splitting terms in Eq. ( 47 1 have explicitly included 



photon radiation by the quarks. These PDFs are 
relevant for the quark distribution prior to a hard 
interaction. Thus, it would be double-counting 
if one included radiative corrections for a quark 
prior to a hard interaction, since these represent 
the same terms that were included in Eq. ( |47| . 
Such a procedure corresponds to the "DIS fac- 
torization" scheme, which assumes that the O(a) 
corrections arising from photon emission from in- 
coming quarks are included in the definition of the 
quark PDFs. The consistent treatment of partonic 
radiative corrections is discussed in some detail by 
Dicncr, Dittmaier and Hollik 87. 



Experimental Limits on Valence Quark 
Charge Symmetry 



There have been no direct observations of any 
violation of partonic charge symmetry. As a re- 
sult we have at present only upper limits on the 
magnitude of parton CSV. We also have the indi- 
rect evidence for partonic CSV from the global fits 
carried out by the MRST group and discussed in 



Sect. Ill A From Eq. (29) in Sec. II D we can ob 



MRST results for valence quark CSV are quanti- 
tatively quite similar to the Gliick analysis. 



tain a relation between the F 2 structure function 
in charged-lepton DIS, and the average of the F 2 
structure functions for neutrinos and antineutri- 
nos, both on isoscalar targets. The difference be- 
tween these two (appropriately normalized) struc- 
ture functions is given by two components. The 
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first is a contribution from strange and charmed 
quarks, and the second is a contribution from par- 
tonic CSV. 

At small Bjorken x, comparison of these two 
structure functions will provide a linear combina- 
tion of heavy quark PDFs and charge symmetry 
violating parton distributions. Extracting limits 
on parton CSV then requires very accurate knowl- 
edge of heavy quark parton distributions. This is 
further complicated by the fact that sea quark dis- 
tributions increase quite rapidly at very small x; 
so the fractional contribution of partonic CSV is 
likely to be small in this region. 

One has two possibilities for placing stronger ex- 
perimental limits on partonic CSV. The first is to 
go to large Bjorken x. Since the heavy quark PDFs 
are quite small at large x, the relative contribution 
of charge symmetry violating parton distributions 
will be significantly larger. The second possibil- 
ity is to take the first moment by integrating the 
parton distributions over all x. One then uses the 
valence quark normalization condition. This con- 



muon and neutrino F2 structure functions, when 
corrected for heavy quark contributions and the 
factor 5/18, should be one independent of x and 
Q 2 , in the naive parton model. The factor 5/18 in 
Eq. (p9| and in Eq. ([52| reflects the fact that the 



virtual photon couples to the squared charge of the 
quarks while the weak interactions couple to the 
weak isospin. The quantity R c is sometimes called 
the "charge ratio," and the relation between the 
F2 structure functions is often termed the "5/ 18 th 
rule." In Eq. ( [52] ) the final line is expanded to 
lowest order in the (presumably small) CSV terms. 

The quantity R c {x) requires knowledge of the 
heavy quark PDFs. For example, the observ- 
ables most sensitive to strange quark distributions 
are cross sections for opposite sign dimuon events 
produced from neutrino DIS on nuclei |621 163] . 
Once the strange quark distributions have been ex- 
tracted from the dimuon production process, they 
can be inserted in Eq. ( 52 1 . The intrinsic charm 



dition, given in Eqs. (17 1 and (18), requires that 



PDFs are generally quite small; however a signifi- 
cant amount of data is collected near charm quark 
threshold, where it is important to take proper 
account of the charm mass. Comparing the F2 
structure functions for lepton-induced processes 
with the F 2 structure functions from weak pro- 
cesses mediated by IF-exchange, one can in princi- 
ple measure both the magnitude and x dependence 
of parton CSV. Clearly, since extraction of par- 
ton CSV distributions depends on precise knowl- 
edge of strange and charm PDFs, our knowledge of 
these quantities will be strongly correlated. Cer- 
tainly this is the case at low x, where the sea quark 
distributions (including strange quarks) are large. 

The charge ratio provides the strongest direct 
limits to date on parton CSV. There should be 
no additional QCD corrections to this relation so 
it should be independent of Q 2 , provided that the 
structure functions are calculated in the so-called 
"DIS scheme," where the F2 structure functions 
sufficiently high energies we can construct the ratio are defined to have the form F 2 (x) = x £\ e 2 qf(x) 

to all orders, where is the quark charge appro- 
priate for either the electromagnetic or weak inter- 
actions. For example, the CTEQ4D parton distri- 
butions [5T] were determined in the DIS scheme. 

Ever since one has been able to extract the F2 
structure functions and hence the parton distribu- 
tion functions from both muon and neutrino DIS, 
one has had the possibility of constructing the 
charge ratio using Eq. ( [52] ) . Within error bars, 
the results have always been consistent with the 
assumption of parton charge symmetry. However, 
until a few years ago the charge ratio gave only 
qualitative upper limits on CSV, because of the 



the first moment of the valence quark CSV terms 
vanish. Consequently, if one integrates parton dis- 
tributions over all x the valence quark CSV terms 
must give zero. Following this integration the only 
remaining contributions will be from the first mo- 
ment of the sea quark CSV. 



1. The "Charge Ratio:" Comparison of F2 Structure 
Functions in Reactions of Muons and Neutrinos 

Currently, the strongest upper limit on parton 
CSV distributions is obtained by comparing the 
i*2 structure functions measured in CC reactions 
induced by v and V, and the F2 structure function 
for charged lepton DIS, both measured on isoscalar 
targets. Using the relation derived in Eq. (29), at 
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In Eq. ( 52 1 the function F 2 (x) is the average of 



the CC F2 structure functions induced by v and 
v and defined in Eq. ( 30 ) . In the denominator 
of the second line of Eq. ( p52"| the sum is taken 
over all quark flavors. Eq. fl52j) shows that in the 
limit of exact charge symmetry, the ratio of the 
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great difficulty in obtaining precise absolute neu- 
trino cross sections, and because of the number 
of corrections that must be taken into account. 
These corrections include: relative normalization 
between lepton and neutrino cross sections; con- 
tributions from strange and charm quarks; higher 
twist effects on PDFs; and heavy quark thresh- 
old effects. In addition, one must be able to 
separate F2 and F$ structure functions in u CC 
reactions. Another potentially important effect 
is heavy target corrections in neutrino reactions. 
The most precise lepton structure functions are 
obtained from deuteron targets, while the most ac- 
curate neutrino cross sections are extracted from 
experiments on heavy targets like iron, so it is nec- 
essary to correct the neutrino F% structure func- 
tions for heavy target effects, and also for effects 
arising from the fact that iron is not an isoscalar 
target. These effects include shadowing and an- 
tishadowing at small x < 0.1, EMC effects for 
0.2 < x < 0.6 EOJ EI] and Fermi motion at large 
x. 

Earlier analyses compared the muon F2 struc- 
ture functions of Meyers et al. on iron |llj to F% ob- 
tained from CCFRR neutrino measurements IT2"1. 



The extracted ratio R c of Eq. ( 52 ) was consistent 



with unity, except possibly at the largest value 
x = 0.65. The experimental data was consistent 
with zero charge symmetry violation and ruled out 
very large violation of parton charge symmetry. 
However, the extracted charge ratio had errors of 
several percent. Because of the factor of 3/10 in 
the last line of Eq. (52 1, an error of 5% in the 



charge ratio would lead to upper limits on par- 
ton CSV at roughly the 15% level. From our dis- 
cussion of the phenomenological and theoretical 
estimates of parton CSV summarized earlier, we 
expect that the CSV contribution to the charge 
ratio will not exceed a few percent at any value 
of x. Consequently, a measurable deviation of the 
charge ratio from unity, at any value of x, would 
be very interesting but experimentally quite chal- 
lenging. 

In recent years we have obtained significantly 
more precise DIS data for both muons and neu- 
trinos. This should allow us to make more strin- 
gent tests of parton charge symmetry. The NMC 
group [891 E0J EI] measured the F 2 structure func- 
tion for muon interactions on deuterium at en- 
ergy = 90 and 280 GeV. The NMC measure- 
ments are more precise than the earlier BCDMS 
muon scattering results on deuterium [H] and car- 
bon [13] , or the SLAC electron scattering results 
[HITS]. The CCFR group [T7] extracted the F 2 



structure function for v and V interactions on iron 
using the Quadrupolc Triplet Beam at Fermilab. 
They also performed a comprehensive comparison 
of their neutrino data with the NMC muon re- 
sults [171 1 



In Fig. 
fl52|) vs. x. 
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we plot the charge ratio 
The solid circles give the 



R c of Eq. 

charge ratio comparing the NMC and CCFR mea- 
surements. The open triangles give the charge ra- 
tio comparing CCFR with BCDMS data, and the 
solid triangles compare CCFR neutrino data with 
the SLAC electron scattering measurements. 

Analysis of the charge ratio as a function of x 
should in principle provide a test of parton charge 
symmetry in both the valence and sea regime. In 
the region 0.1 < x < 0.4, where valence quarks 
should dominate, the charge ratio is consistent 
with unity, with errors on the charge ratio at about 
the 3% level. From Eq. (52 1, this would provide 



upper limits on valence quark CSV of about 10%. 
For larger values of x the upper limit on errors 
in the charge ratio is in the 5-10% level, due pri- 
marily to the poorer statistics and the large Fermi 
motion corrections that become important at very 
large x. Particularly after including heavy tar- 
get corrections for the ^-iron measurements, the 
charge ratio appeared to deviate significantly from 
one at the smallest values x < 0.1 92J. The de- 
viation appeared to grow with decreasing x and 
reached values as large as 15-20%. This is ap- 



parent from Fig 11 which shows that W c is def- 
initely less than one for small x < 0.1. Boros, 
Londergan and Thomas [20 HU] studied the origin 
of this discrepancy. They suggested that, assum- 
ing that the identification of the neutrino F2 CC 
structure functions was reliable, the most likely 
explanation for this anomaly was a substantial vi- 
olation of charge symmetry in the nucleon sea. 

Boros et al. showed that one needed sea quark 
CSV of at least 25% in order to explain this dis- 
crepancy in the charge ratio. This apparent vio- 
lation of charge symmetry was extremely surpris- 
ing, as it was at least an order of magnitude larger 
than theoretical estimates. Bodek et al. E2J ques- 
tioned whether such a large CSV effect was consis- 
tent with other experiments. They analyzed the 
W boson charge asymmetry obtained in pp exper- 
iments from the CDF group at the Fermilab Teva- 
tron 94J. Since this experiment involves proton- 
antiproton scattering, CSV effects do not enter di- 
rectly. However, Bodek and collaborators argued 
that the most precise PDFs arise from charged lep- 
ton DIS on isoscalar targets. Because the F2 struc- 
ture functions are weighted by the squared charge 
of the quarks, they are most sensitive to the up 
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11: Charge ratio R u c {x) of Eq.[52 
CCFR v 



vs. x. Solid cir- 



Fe data, Ref . [T7j and [i + D measure- 
ments from NMC, Ref. [§3 [90] [91]. Open triangles: 
CCFR ^ data and fi+ D measurements from BCDMS, 
Ref. [13]. Solid triangles: CCFR and SLAC electron 
scattering data, Ref. [101 115] . 



quarks in the proton and neutron. Thus to a sig- 
nificant degree our identification of d p is obtained 
from u n plus the assumption of parton charge sym- 
metry. Bodek examined two different methods for 
extracting CSV distributions from the data, and 
calculated the effect on the W charge asymmetry 
which would arise from CSV effects of the magni- 
tude assumed by Boros et aZ.The results are shown 
in Fig. [12] 
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FIG. 12: [color online] The W charge asymmetry from 
pp reactions at the Tevatron. The experimental points 
are those of the CDF group, Ref. [93] • The solid and 
dash-dot curves represent two fits using the CTEQ4M 
PDFs [81] with no CSV terms. The dashed and dotted 
curves represent two different assumptions by Bodek 
et al, Ref. [93] using sea quark CSV distributions cal- 
culated by Boros et al, Ref. [H[2D]. 



In Fig. [12] the solid and dash-dot curves are fits 
to the CDF data [Hj using the CTEQ4M parton 
distributions [81] with no parton CSV terms. The 
dashed and dotted curves resulted from two differ- 
ent assumptions by Bodek et al. for the large sea 
quark CSV terms of Boros et al. QI1[20]. The CDF 
measurements are very sensitive to the sea quark 
distributions, and Bodek argued that the large sea 
quark CSV was incompatible with those experi- 
mental results. Although Bodek and collaborators 
examined only two potential ways of defining par- 
ton CSV distributions, it is hard to imagine that 
sea quark CSV of this magnitude could be made 
consistent with the W charge asymmetry data. 

This issue was eventually resolved when the 
CCFR collaboration re-analyzed its neutrino data 
[22 and the low- a; discrepancy disappeared. There 
were two primary reasons for this change. The first 
was an improved treatment of charm mass correc- 
tions. This was particularly important for the low- 
x data, which were taken in a region close to charm 
threshold. In analyzing this data it is necessary 
to take into account accurately the charm quark 
mass. The initial analysis used a "slow rescal- 
ing" hypothesis due to Georgi and Politzer [95] [96] 
to account for charm mass corrections. The re- 
analysis involved NLO calculations, which account 
for massive charm production using variable-flavor 
techniques [2T] |97] |98] . 

The second significant effect involved the sep- 
aration of structure functions in charged-current 
v DIS. The sum of v and V charged-current DIS 
cross sections gives a linear combination of Fi and 
F 3 structure functions, 

2(l-y-y 2 /2)F™ N °(x,Q 2 ) 



dxdy dxdy 



(y~y 2 /2)AxF 3 (x,Q 2 ) . (53) 



In Eq. ( 53 1 , the quantity AxF 3 (x) is the difference 



in the F 3 charged-current structure functions for 
neutrino and antincutrino beams. 



AxF 3 (x) = xF 3 w (x) 



xF 3 w (x) ; 



AxF^ix) 



x 2(s + (x) - c + (x)) 



+ Sd+(x)-6u + (x) 



(54) 



The second line in Eq. (54 1 is valid to leading order 



in QCD for an isoscalar target and for sufficiently 
high Q 2 . In these limits, assuming the validity of 
charge symmetry, the quantity AxF 3 is sensitive 
only to heavy quark distributions. 

For simplicity in Eq. ( 53 1 we have dropped terms 



of order M / s and have set the longitudinal to 
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transverse ratio R to zero (these approximations 
were not made in re-analyzing the data). In the 
initial analysis |17j . data for a given x bin was av- 
eraged over all y, and the AxF 3 structure function 
was estimated using phenomenological PDFs. In 
the re-analysis the data was binned in x and y so 
that both F2 and AxF^ could be extracted [22] . 
The experimental values for AxF$ differed sub- 
stantially from the phenomenological predictions. 
From Eq. (53), a change in AxF^ will affect the 



values extracted for the charged-current F 2 neu- 
trino structure functions. The combined effect of 
the NLO treatment of charm production, and the 
model-independent extraction of AxF^ removed 
the small- x discrepancy. The charge ratio R c of 
Eq. [52] is now unity to within experimental error, 
even at small x. 



the magnitude of CSV effects in the sea, for values 
x > 0.015. To obtain more quantitative limits on 
CSV, it will be necessary to obtain reliable esti- 
mates for the few remaining uncertainties in this 
comparison. Perhaps the largest undetermined 
correction remains the shadowing of parton distri- 
butions for v—Fe interactions. Experimental anal- 
yses have assumed that the nuclear shadowing cor- 
rections are the same for neutrinos (virtual W's) 
as for charged leptons (virtual photons). Boros et 
aL[9"2"] showed that one could expect substantially 
different shadowing for W's than for photons, pri- 
marily because the Ws couple to axial currents 
as well as to vector currents. This was further 
expanded by Brodsky et al. \W0\ 1101) who calcu- 
lated both shadowing and anti-shadowing effects 
for neutrino DIS. 
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FIG. 13: The ratio 5F%/(18F£), calculated by the 
CCFR Collaboration, Ref. [22j . Curves are for vari- 
ous NLO parton calculations. Solid curve: Ref. |97| : 
dotted curve: Ref. [98]; dash-dotted curve: Ref. |99j . 
Neutrino structure functions from the CCFR group, 
Ref. [22]. Solid points, NMC muon data, Ref. [91]; 
open circles: BCDMS data, Ref. [14] : diamonds: 
SLAC data, Ref. Q3]. 



The results are shown in Fig. [13] These graphs 
plot the ratio 5F£/(18F£) from the CCFR re- 
analysis, vs. Q 2 for various values of x. The dif- 
ferent data points involve muon DIS experiments 
from NMC, BCDMS, and SLAC [H HI [T5] EE] . 
The curves are NLO analyses using various meth- 
ods for including charm mass effects [2T] [97] l98l 
199] . The previous low- a; discrepancy between the- 
ory and experiment has largely disappeared. This 
allows one to place qualitative limits of < 10% on 



2. Charge Symmetry and Determination of the 
Weinberg Angle 

The NuTeV group [25] US have measured total 
charged-current and neutral-current cross sections 
for v and v on an iron target. From these measure- 
ments they made an independent determination of 
the Weinberg angle, motivated by a procedure ini- 
tially suggested by Paschos and Wolfenstein [102] . 
Paschos and Wolfenstein showed that a ratio of 
total cross sections for neutral current (NC) and 
charged current (CC) interactions for v and V on 
an isoscalar target Nq gave the remarkably simple 
Paschos- Wolfenstein (PW) relation 



R- 



uN 

JVC 



Pi 



1 • 2 

= sin 

2 



(55) 

In Eq. (55 1, the quantities are the total NC and 
CC cross sections for v and V on an isoscalar tar- 
get, and the quantity po = Mw/(Mz cos 9w) is 
one in the Standard Model. The brackets denote 
integration of the cross sections over all Bjorken 
x. Although the individual cross sections depend 
upon details of parton distributions, the ratio of 
these combinations contains no dependence upon 
PDFs, and in addition a number of experimental 
effects cancel. 

The NuTeV collaboration used the Sign Selected 
Quadrupolc Train beamline at Fermilab to sep- 
arate v and V arising from pion and kaon de- 
cays following the interaction of 800 GeV pro- 
tons. The resulting interaction events were ob- 
served in the NuTeV detector, where they were 
required to deposit between 20 GeV and 180 GeV 
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in the calorimeter. CC and NC events were dis- 
tinguished by the event length in the counters, as 
CC events contained a final muon that penetrated 
substantially farther than the hadron shower. The 
NuTeV collaboration measured the individual ra- 
tios R v and R v defined by 



r uN 



R 



R 



(56) 



In terms of the ratios defined in Eq. ( 56 1 , the PW 
ratio has the form 



R~ 



R v -rR v 
1 - r 



(57) 



The NuTeV group measured R v = 0.3916 ± 
0.0007 and R v = 0.4050 ± 0.0016. The quantity 
r = 0.499±0.005 was taken from the world average 
of v — Fe charged-current DIS reactions |1031 1104j 
and from measurements by the CCFR collabora- 
tion [17]. Since acceptances and cuts differ for 
v and V reactions, they did not directly construct 



the Paschos- Wolfenstein ratio via Eq. ( 55 1 ; instead 



the measured NC/CC ratios were compared with 
a Monte Carlo simulation of the experiment, from 
which they extracted the on-shell value for the 
Weinberg angle sin 2 6 W = 0.2277± 0.0013 (stat) ± 
0.0009 (syst). This value is three standard devia- 
tions above the measured fit to other electroweak 
processes, sin 2 6 W = 0.2227 ± 0.00037 [3T]. 

In a given renormalization scheme, the effec- 
tive Weinberg angle sin W will acquire n Q - 
dependence from radiative and loop corrections 
[TM1 [TfJrJl HOT] . In Fig. [I4]we plot the effective 
value for sin 2 9 W vs. Q with the results of sev- 
eral experiments. The curve is that of Erler and 
Ramsey-Musolf |105j . calculated in MS scheme. 
The experimental points represent atomic parity 
violation in Cesium (APV) [108 , a M0ller scatter- 
ing measurement from experiment El 58 at SLAC 
(Qw(e)) [105], a series of measurements at the Z 
pole at LEP and SLD (Z-pole) [21], the forward- 
backward asymmetry of e + e~ pairs produced from 
pp collisions at the Tevatron by the CDF group 
(A FB ) [TT0], and the NuTeV result (i/-DIS) [U]. 
Additional points in red give error estimates for 
two proposed experiments, e — D parity- violation 
in DIS (eD-DIS) jl 1 lj . and an approved Q-weak 
experiment involving DIS for polarized electrons 
on protons at Jefferson Laboratory (Q w (p)) |112j . 
Although the error bars on some of the experi- 
ments are rather large, the E158 experiment estab- 
lishes the Q 2 dependence of the effective Weinberg 



angle at the 6a level, and the NuTeV experiment 
is claimed to differ from the expected value by 3a. 
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FIG. 14: [color online] Effective value for sin u\y vs. 
Q. Curve is from Ref. [105] . calculated in MS scheme. 
Experimental points represent atomic parity violation 
(APV) [108] , M0ller scattering (Q w (e)) [109] . measure- 
ments at the Z pole [3T], NuTeV (z/-DIS) [H], and 
forward-backward asymmetry from CDF (Afb) |1 10J ■ 
Points in red represent error estimates for two pro- 
posed experiments, e — D parity violating scattering 
(eD-DIS) [111] , and the Q-weak experiment (Q w (p)) 



The NuTeV result, which implies an effective 
left-handed coupling of light quarks to the neu- 
tral current that is about 1.2% smaller than ob- 
tained from other electroweak data, is rather sur- 
prising. The status of what has been termed the 
"NuTeV anomaly" has recently been summarized 
|113j . Davidson et al. |114j considered a number 
of corrections from physics outside the Standard 
Model. It is quite difficult for new physics to ex- 
plain the NuTeV finding since such effects have to 
agree with both the NuTeV result and also with 
the very precise measurements of EW effects at 
LEP, which constrain some parameters to a few 
parts per thousand. As a result, most recent ef- 
forts have focused on effects within the Standard 
Model. At present the three most likely "QCD ef- 
fects" are: effects due to radiative corrections or 
nuclear effects in the neutrino reactions; contribu- 
tions from strange quark momentum asymmetry; 
or charge symmetry violation in parton distribu- 
tions. 

To lowest order in the strong coupling a s , one 
can calculate various analytic corrections to the 
Paschos-Wolfenstein relation. These can be writ- 
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ten in the form 



5R- = 



(-1 + |sin 2 W ) 

U v + D v 
g _ SD V -SU V 



(N-Z) 
A 



(U v - D v ) 



(58) 



Eq. ( 58 ) gives estimates of the corrections to the 
NuTeV result. Terms of the form Q v denote the 
second moment (integral over all x) of a given fla- 
vor valence distribution, e.g., U v = (x[u(x)—u(x)]) 
represents the total fraction of the proton momen- 
tum carried by up valence quarks. The first term 
is an isoscalar correction due to excess neutrons 
in the iron target. The NuTeV collaboration have 
taken this correction into account. The correction 
is large (of order -.008), but should be known to 
within a couple percent. The second and third 
terms represent respectively contributions from a 
possible strange quark momentum asymmetry and 
from parton charge symmetry violation. 

Because the NuTeV group did not directly con- 
struct the Paschos-Wolfcnstcin ratio, Eq. d58l 



gives only an estimate of the effects of these con- 
tributions to the NuTeV experiment. The NuTeV 
group has provided functionals that give the sensi- 
tivity of their experiment (in Bjorken x) to various 
quantities |30| . e.g., charge symmetry violation or 
a strange quark momentum asymmetry. To ob- 
tain a quantitative result for a particular effect, 
one multiplies the effect in question by the ap- 
propriate functional and integrates over x. For 
example, corrections to the Paschos-Wolfenstein 
relation depend only on valence quark properties; 
sea quarks give no contribution to that relation. 
Sea quark corrections to the NuTeV experiment 
are much smaller than the corresponding valence 
quark contributions, but they are not zero. 

Radiative corrections, which involve coupling of 
soft photons to the final muon line, are impor- 
tant for CC events, and constitute a substantial 
correction. Recently, Dicner, Dittmaier and Hol- 
lik have re-calculated the radiative corrections, in- 
cluding all corrections of order O(a), and a num- 
ber of additional higher-order corrections [57J HE] ■ 
They find some differences from the older radia- 
tive correction program of Bardin and Dokuchaeva 
[115] . Diener et al. also include new terms that 
result from including electromagnetic coupling in 
the QCD evolution equations. Such terms have re- 
cently been included by the MRST group [24], and 
by Gluck et al. |25j : these contributions were re- 
viewed in Sect. IIII B 11 Diener et al. estimate that 
radiative correction effects would remove about 



one-fourth of the NuTeV anomaly [BE]- Note 
that these corrections are renormalization scheme 
dependent. The NuTeV group is currently re- 
analyzing their data, using a radiative corrections 
code provided by Diener et al. [116J . 

The NuTeV measurements require nuclear cor- 
rections for the structure functions. Kumano 
[55] calculated a "modified PW relation" for nu- 
clei, and Hirai et al. |55j estimate that nuclear ef- 
fects could remove up to one-third of the NuTeV 
anomaly. They assume that nuclear shadowing for 
neutrinos is identical to that for charged leptons 
[53]. Kulagin and Petti [551 157] also considered 
nuclear effects, particularly in neutrino deep in- 
elastic scattering reactions. Miller and Thomas 
[117] point out that v shadowing effects could dif- 
fer significantly from shadowing of muons |92] , 
They also emphasize that shadowing produces dif- 
ferent effects for CC and NC events. Brodsky et 
al. |100l 1101] made a detailed calculation of shad- 
owing and antishadowing arising from multigluon 
exchange. They conclude that nuclear shadow- 
ing effects could account for roughly 20% of the 
NuTeV anomaly. 

Another effect that might contribute to the 
NuTeV result arises from a strange quark mo- 
mentum asymmetry; this is the second term in 
Eq. (58 1. As was discussed in Sect. II B it is pos- 



sible that s(x) ^ s(x). Although the first moment 
of s — s must be zero (there is zero net strangeness 
in the proton), the second moment 



S~ 



(x[s(x) — s(x)]) = (xs (x)) 



(59) 



(see Eqs. ( 14 1- ( 16 1 ) need not vanish. A nonzero 
value for S~ would mean that the net momentum 
carried by strange quarks and antiquarks was un- 
equal. From Eq. (58 1 we see that if the strange 
quark momentum asymmetry S~ is positive (in 
this case, strange quarks would carry more of the 
nucleon's momentum than strange antiquarks), 
this would decrease the extracted value of sin 2 9yy , 
and decrease the discrepancy with the expected 
value of the Weinberg angle; conversely, a nega- 
tive value of S~ would increase the discrepancy. 

The most direct knowledge of strange quark dis- 
tributions comes from measurements of opposite 
sign dimuons produced in neutrino-induced nu- 
clear reactions. In such reactions, dimuon produc- 
tion from v (v) beams is sensitive to the s (s) dis- 
tribution, so that in principle comparison of these 
cross sections could enable one to determine differ- 
ences between s and s PDFs. These cross sections 
have been extracted by the CCFR [55] and NuTeV 
[55] collaborations. In the CCFR experiment the 
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FIG. 15: [color online] The quantity xs (x) — x[s(x) — 
s(x)], vs. x, as extracted by the NuTeV Collaboration, 
Ref. [BS]. Values are obtained for Q 2 = 16 GeV 2 . The 
outer error band is the combined error, while the inner 
band is without the uncertainty in the semileptonic 
branching ratio B c . 



v and V beams are not separated and the type of 
reaction is inferred from the charge of the faster 
muon, while the NuTeV experiment uses separated 
v and V beams. 

For some time there was disagreement as to the 
interpretation of the dimuon experiments and ex- 
traction of the strange quark PDFs. The NuTeV 
group analyzed the dimuon cross sections and ex- 
tracted strange distributions [3U]. Their results 
were consistent with a small value for s~(x), with 
a second moment that was zero or slightly negative 
[50] ; the value that they extracted would increase 
the discrepancy in the Weinberg angle to about 
3.7a. On the other hand, the CTEQ group [64] 
estimated that S~ was most likely positive, and 
they suggested that this could remove roughly 1/3 
of the NuTeV anomaly. The CTEQ global analy- 
sis of s~(x) was dominated by the CCFR/NuTeV 
data [621 163j for opposite-sign dimuon production 
in neutrino DIS, so it was unclear why the two 
groups obtained differing results. Since then the 
CTEQ and NuTeV groups have collaborated on 
the data analysis, recently obtaining consistent re- 
sults. 

The latest NuTeV result obtained by Mason et 
al.\EE\ yields a best value for S~ that is positive. 



Fig. 15 plots the quantity xs (x) vs. x from the 



latest NuTeV analysis. This results in a quantity 



In Eq. (60), the quantity "external" refers to 



the contribution due to uncertainties on external 
measurements. The strange quark asymmetry of 
Eq. ( 60 ) would remove roughly one-third of the 



NuTeV anomaly. The NuTeV group provides a de- 
tailed error analysis of the quantity S~ . It is quite 
sensitive to two quantities. First is the semilep- 
tonic branching ratio B c ; the outer band in Fig. |15| 
shows the result for S~ with the B c uncertainty, 
and the inner band is the result without the B c 
uncertainty. The second is the point at which the 
quantity xs~(x) crosses zero (it must cross zero at 
least once so the first moment of s — s is zero). The 
current best fit crosses zero at a very small value 
x ~ 0.004. This means that the quantity s~(x) 
would have a very large negative spike at very low 
x. It is very difficult to imagine a physical mech- 
anism that would cause s~(x) to change sign at 
such a small value of x. 

If one allows the zero-crossing point to increase, 
then the resulting value of S~ decreases, but the 
X 2 value also increases somewhat. The best value 
obtained by Mason et al. S~ — 0.00196 occurs for 
a zero crossing of x = 0.004 and \ 2 = 38.2 for 
37.8 effective degrees of freedom. Fig. 



16 



gives ex- 
amples of the relation between the zero-crossing 
point, the resulting curve of s~{x) vs. x, and the 
resulting \ 2 - For example, when the zero cross- 
ing moves to x = 0.15, then one obta ins S~ = 
0.00007 but x 2 increases to 53.4. Fig. Ilj] shows 



the strong correlation between strange quark mo- 
mentum asymmetry s~ (x) and the zero-crossing 
point. 

The contribution from charge symmetry violat- 
ing parton distributions to the NuTeV anomaly 
(the last term in Eq. (58)) can be estimated by 



folding quark CSV distributions with the function- 
als provided by the NuTeV group. This contribu- 
tion is dominated by valence CSV distributions. 
Using the phenomenological CSV PDFs obtained 
by the MRST global fit [55], valence CSV with 
k = —0.6 would completely remove the NuTeV 
anomaly, whereas the value k = +0.6 would make 
it twice as large. Both of these values are within 
the 90% confidence level in the MRST global fit. 
Thus the uncertainty in parton charge symme- 
try violation as calculated by MRST is capable of 
removing completely the NuTeV anomaly in the 
Weinberg angle. 

We can also investigate parton CSV contribu- 
tions to the NuTeV result from theoretical calcu- 



S~ = 0.00196 ± 0.00046(siai) ± 0.00045(sysi) rations. From Eq. (58), the contribution from CSV 
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FIG. 16: [color online] The quantity xs~ (x) — x[s(x) — 
s(x)], vs. x, as extracted by the NuTeV Collaboration, 
Ref. [5S]. Three different results are shown, corre- 
sponding to different values of the zero-crossing point. 
The x 2 value is listed for each curve. 



to the PW ratio is proportional to 



6R- 



SD V - SU V 
U v + D v 



(61) 



Thus the contribution from CSV to the PW ratio 
is related to the second moment of the CSV va- 
lence PDFs, divided by the total momentum car- 



ried by up and down valence quarks. In Sect. Ill B 
we showed that Sather's analytic approximation 
for valence charge symmetry violation gave an an- 
alytic expression for the second moment of these 
distributions (see Eq. ( 44 1 ) . Since these quanti- 



ties depend only on total momentum carried by 
valence up and down quarks, quantities which 
are reasonably well determined, it was argued 
that the second moments of valence parton CSV 
were essentially model-indepedent quantities [86 . 
This partonic CSV correction would decrease the 
anomaly in the PW ratio by roughly 40%. 

However, as we have stated the NuTeV group 
did not measure the Paschos-Wolfenstcin ratio. If 
instead one uses the theoretical CSV distributions 
from Rodionov et al. |74j with the functionals pro- 
vided by NuTeV, then one finds that valence par- 
ton CSV removes about 1/3 of the anomaly in 
sin 2 6\y [77J [SS] . Charge symmetry violation aris- 
ing from the QED splitting mechanism described 
in Sect. IIIB1 would remove another 1/3 of the 
anomaly. 

Davidson and Burkardt |118j have estimated the 
effect on the Paschos-Wolfenstein relation arising 
from nuclear charge symmetry violation, i.e. the 
fact that protons are more weakly bound than 



neutrons due to Coulomb effects. Their results 
suggest that these nuclear Coulomb effects would 
increase the magnitude of the NuTeV anomaly by 
roughly 20%. A new nuclear mechanism has re- 
cently been suggested by Cloet et al. [119] . We 
call this "pseudo CSV," since it produces effects 
that mimic those arising from charge symmetry 
violation. We discuss this effect and the implica- 
tions for the NuTeV experiment in the following 
section. 

We have shown that it is necessary to consider 
a number of "QCD effects" within the Standard 
Model, in order to obtain precise results for the 
NuTeV experiment. Small but non- negligible con- 
tributions are likely from nuclear effects on par- 
ton distributions and strange quark effects. At 
this time, charge symmetry violation appears to be 
the only mechanism capable of single-handedly re- 
moving the entire NuTeV anomaly. An additional 
nuclear effect which we term "pseudo CSV" is dis- 
cussed in the following section; this effect is ca- 
pable of making a substantial contribution to the 
NuTeV measurement. Another possibility would 
be that a new treatment of radiative corrections 
might produce significant corrections to the ex- 
tracted value for the Weinberg angle. However, a 
re-analysis of the NuTeV data using the newer ra- 
diative corrections [57J EE] has not been published 
at this time |116j . 



D. Pseudo-CSV Nuclear Effects 

Many of the tests of partonic charge symme- 
try violation, and some applications that rely on 
CSV have been carried out with neutrino beams 
and often with an Fe target - simply to increase 
the event rate. The NuTeV anomaly, as an ex- 
ample, was carried out with an Fe target, even 
though the Paschos-Wolfenstein relation is only 
valid for an isoscalar target. Of course, the cross 
sections for v and v scattering were corrected for 
the small number of excess neutrons. However, 
as pointed out in a recent paper by Cloet, Bcntz 
and Thomas |119j . this will in general not be suf- 
ficient. It has been understood for some time [55] 
that the famous "EMC effect" [FUE], the nuclear 
modification of the Fi structure function in elec- 
tromagnetic DIS reactions, cannot be understood 
simply in terms of the Fermi motion and binding of 
free nuclcons, but the actual quark structure of the 
bound nucleon must also be modified in a signifi- 
cant way. A number of relatively successful models 
have been constructed |120l 11211 1122] , based upon 
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the self-consistent modification of the bound nu- 
clcon structure in the relativistic mean scalar and 
vector potentials generated in a nuclear medium. 
The new realization in the case of Fe and, indeed 
any other nucleus with N ^ Z, is that there will be 
an isovector piece of the EMC modification of the 
bound nucleon structure associated with the extra 
neutrons. Most important, this effect will modify 
the structure of all of the neutrons and protons in 
the nucleus, not just the excess neutrons. 

As the dominant piece of the isovector interac- 
tion in a relativistic mean field theory is usually as- 
sociated with the p meson, it will have a Lorentz 
vector character, with the (i-quarks feeling more 
repulsion and the M-quarks more attraction. For 
this reason the sign of the effect is exactly the same 
as that found in the calculations of CSV which we 
have described earlier. If one ignores this medium 
modification it will appear as though the CSV is 
enhanced in a nucleus with N > Z. We stress 
that there is no violation of charge symmetry - 
the isovector interaction is completely consistent 
with isospin invariance - but to an observer un- 
aware of the isovector EMC effect it will appear 
like CSV. An estimate of the impact of this ad- 
ditional EMC effect on the NuTeV analysis [TT9] 
based on a nuclear matter calculation, reduces the 
NuTeV result for sin 2 9 W from 0.2277 to 0.2245, 
within 1 a of the Standard Model value. If the 
effect of true CSV is added [77J |53] (as was dis- 
cussed in the preceding section), the discrepancy 
is removed entirely. 
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FIG. 17: [color online] The EMC ratio {x) / (x) 
vs. x, at a value Q 2 — 10 GeV 2 , predicted by the 
model of Cloet, Bentz and Thomas |119] as a function 
of the proton/neutron ratio Z/N < 1. Solid curve: 
Z/N — 1; dashed curve: Z/N = 0.8; dash-dot curve: 
Z/N = 0.6; triple dot-dashed curve: Z/N = 0.2; dot- 
ted curve: neutron matter (Z/N — 0). 

It will clearly be very important to look for spe- 



cific processes which could confirm this theoretical 
analysis of the isovector EMC effect. This model 
predicts a significant and characteristic A depen- 
dence of the ratio of the nuclear F 2 electromag- 
neti c str ucture function with that for the deuteron. 
Fig. [lTjshows the EMC ratio, F£(x)/Fg{x), vs. x 
at 0* = 10 GeV 2 , for various values of N > Z. For 
a neutron excess, the medium modification of the 
u quarks should be enhanced by coupling to the p° 
field, while the d quark distribution should be less 
modified. For small neutron excess the EMC ef- 
fect, which is initially dominated by the u quarks, 
increases. However eventually the d quark distri- 
bution dominates and the EMC ratio is predicted 
to decrease in the valence quark region. For ex- 
ample, in Au where N ~ 1.5 Z, a very large differ- 
ence is predicted between the ratio of u(x) in Au 
to that in the deuteron, compared with the same 
ratio for (i-quarks. This could be investigated in 
experiments at Jefferson Laboratory following the 
12 GeV upgrade. 

If the "pseudo-CSV" nuclear effect outlined 
here is confirmed experimentally, then it would 
seem that rather than presenting evidence for new 
physics beyond the Standard Model, the NuTeV 
result rather confirms in a fairly dramatic fashion 
the concept that the partonic structure of a bound 
nucleon is modified in a profound way. 



E. Dedicated Experiments Sensitive to 
Valence Quark Charge Symmetry 

In the preceding Section, we reviewed existing 
experiments and showed the limits they placed on 
charge symmetry and flavor symmetry violation 
in parton distributions. In this section we discuss 
various dedicated experiments that might tighten 
the limits on parton charge symmetry, and we re- 
view the conditions that would be necessary in 
order that these experiments could detect parton 
CSV at levels that are allowed from current phe- 
nomenological limits. 

As will become clear, the experiments described 
are looking for quite small effects, of the order of 
one or a few percent. In addition, one generally 
has additional terms arising from other effects such 
as heavy quark contributions. These must be un- 
der control before one can isolate parton CSV ef- 
fects. Finally, several of these require subtracting 
cross sections from two separate measurements. 
These experiments are then very sensitive to rela- 
tive normalizations. Although these are not easy 
experiments, it is also true that even tighter up- 
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per limits on CSV contributions to parton distri- 
bution functions could improve dramatically our 
understanding of these effects. 



1. Drell-Yan Processes Initiated by Charged Pions 

A suitable probe for charge symmetry effects 
should differentiate between up quarks in the pro- 
ton and down quarks in the neutron. This can be 
accomplished by comparing Drell-Yan (DY) pro- 
cesses induced by charged pions on isoscalar tar- 
gets. Drell-Yan processes |123|, I124j proceed via a 
quark (antiquark) from the projectile annihilating 
an antiquark (quark) of the same flavor from the 
target, producing a virtual photon that eventually 
decays into a pair of oppositely charged muons 
with large Q 2 . The process is shown schemati- 
cally in Fig. [18] The left figure shows the pro- 
cess for NN Drell-Yan processes. The right figure 
shows a schematic diagram for ir + — p Drell-Yan 
processes, in the kinematic regime where valence 
quarks dominate for both the pion and nucleon. 
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FIG. 18: Schematic picture of the Drell-Yan (DY) pro- 
cess, production of a fi + — fiT pair with high invariant 
mass through a virtual photon. Left: NN DY process; 
a quark in one nucleon annihilates with an antiquark 
of the same flavor in the second nucleon. Right: % + p 
DY process, in the valence-dominated region of x for 
both 7r + and p. 

We review here the calculation of Londergan et 
al. |1251 1126j , who suggested using Drell-Yan pro- 
cesses initiated by pions to study partonic CSV. 
At large momentum fraction x, the nucleon distri- 
bution is dominated by its three valence quarks, 
while at large x^ the pion is predominantly a va- 
lence q — q pair. For DY processes induced by 
charged pion beams on nucleon targets, in the 
kinematic region of reasonably large Bjorken x for 



both projectile and target quarks, the annihilating 
quarks will come predominantly from the nucleon 
and the antiquarks from the pion. The tt + con- 
tains a valence d (and will annihilate a d quark in 
the nucleon) and n~ a valence u (and will annihi- 
late a nucleon u quark) . 

Therefore, comparison of ir + and 7r~ induced 
DY processes on an isoscalar target such as the 
deuteron should provide a sensitive method for 
comparing d and u valence distributions in the nu- 
cleon. As the x and x^ values of interest for the 
proposed measurements are large, a beam of 50 
GeV pions will produce sufficiently massive dilep- 
ton pairs that the Drell-Yan mechanism is appli- 
cable. A flux of more than 10 9 pions/second is 
desirable. These experiments might be feasible 
for fixed target experiments using the Fermilab 
Main Injector |127j . Alternatively, such exper- 
iments would be possible in the COMPASS ex- 
periment at CERN |128j . provided that one used 
charged pion rather than muon beams. There ex- 
ist some data for 7r + Drell-Yan scattering nuclear 
targets dating from about thirty years ago, Fermi- 




lab experiment E444 [129] and CERN experiment 
WA39 |130j . As a general rule Drell-Yan experi- 
ments with 7r + beams are more difficult than 7r~, 
since the pions arc generally secondary beams aris- 
ing from proton bombardment and one must be 
able to separate the ir + from protons. In addi- 
tion, the DY cross sections for ir~ will generally 
be larger than the corresponding DY cross section 
induced by 7r + , as seen from Eq. ( 62 ) . 



Consider the DY process for a charged pion on 
a deuteron target. Neglecting for the moment sea 
quark effects, at sufficiently large x and i, the 
7r ± -D DY cross sections will be proportional to: 



1 , *, 



a% D (x,x v ) ~ -(dP(x) + d n (x))d (x w ) , 



a^I D (x,x r ) ~ -(u*(x)+u n (x))W (a*). 



(62) 



Consider the ratio R^(x,x n ), defined by 



t>DY 
K ttD 



-DY 



(x, Xtt) 



DY 
T tv+D 



(x, Xyr) 



This ratio will be sensitive to charge symmetry vi- 
olating (CSV) terms in the nucleon valence parton 
distributions. Since theoretical CSV effects are no 
greater than a few percent, sea quark contribu- 
tions for both nucleon and pion must be included. 
To first order in small quantities the DY ratio for 
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pions can be written |1251 1131] : 

, 27T S (X W ) 



Rcs(x) 



n v (x n ) 

4(Sd v (x) - <5u v (ie)) 
3(it v (x) + d v (x)) 



(64) 



The term R sv (x,x 7T ) in Eq. (|64j contains sea- 
valence interference terms which are given in 
Ref. |m 

Eq. (64 1 used charge conjugation invariance and 



assumed charge symmetry for the pion PDFs. In 
this case one can write 



ir v (x) 
tt s (x) 



uf(x) = df{x) = dl (x) = < (x) 
qf (x) = qf (x) = qf (x) = qf (x) 



[q = u, d] 



(65) 



Eq. (64) is valid at sufficiently large x and %„. It 



is expanded to lowest order in both sea quark and 
CSV terms. In Eq. (64 1 a pion CSV term has been 



is not sensitive to differences between parton dis- 
tributions in the free nucleon and in the deuteron 
[T32"l EM CE31 E3S], provided that both neutron 
and proton parton distributions are modified in 
the same way . 

DY i 



In Fig. [P9| we show the ratio R^(x,x v ) vs. x 
for x n — 0.4. For the nucleon PDFs, we used 
the MRST global fit distributions including a va- 
lence CSV contribution of the form of Eq. (35 1. 



The pion PDFs were taken from those of Sutton 
et al. [133]. These were fit to older pion DY NA10 
and E615 experiments [137, 138J ; these pion distri- 
butions can be evolved to higher Q 2 using inter- 
polating matrices supplied by the MRST group. 
Ref. 136 provided several different pion PDFs; for 
Fig. 19 we used pion PDFs for which 10% of the 



pion momentum was carried by the sea. The pion 
and nucleon PDFs were evolved to a typical value 
Q 2 = 25 GeV 2 . In Fig. 19 the solid curve corre- 



sponds to zero CSV contribution, the dashed curve 
is for k = +0.65, and the dotted curve to k = —0.8. 
These curves represent the 90% confidence limit on 
the CSV distributions for the MRST global fit to 




FIG. 19: Theoretical estimate of nucleon CSV term 
Res of Eq. |64l, vs. x, for x n = 0.4 and Q 2 = 25 GeV 2 
from Ref. [131| . Solid curve: no CSV terms, k = 0; 
dashed curve: k — +0.65; dotted curve: k — —0.8. 



given in Eq. (351. The 



neglected; theoretical models predict a very small 
effect from this term |125j . 

The relative DY fluxes for charged pions can 
be obtained by measuring the yield of J/f^'s from 
7r + — D and tt~ — D, which should be identical to 
within one or two percent. The nucleon CSV term 
Rcs(x) in Eq. (64 1 is a function only of x. A num- 
ber of systematic errors should cancel in taking 
the ratio of cross sections. In particular, Eq. (64 1 m Sect. |III A 



valence quark CSV 
CSV contribution is surprisingly large. At x — 0.5 
the limit of the two CSV terms represents about 
a 50% correction to the ratio, while at x — 0.8 the 
contribution is nearly 100%. 

This very large contribution from charge sym- 
metry violation is almost certainly an artifact of 
the fact that the MRST CSV PDFs are indepen- 
dent of Q 2 , while the parton distributions in the 
denomin ator d epend upon Q 2 . This was discussed 
At large values of Q 2 , 



DGLAP 

evolution causes valence parton distributions to 
move to progressively smaller x values. For val- 
ues x > 0.3, the numerator (Q 2 independent) will 
remain large while the denominator will become 
progressively smaller. We expect that the ratios 
shown in Fig. [19] would decrease substantially if 
the CSV parton distributions were evolved in Q 2 . 

In a Drell-Yan ir — D experiment, one would first 
measure DY cross sections over a wide kinematic 
region, and extract the pion valence and sea distri- 
butions. One would then construct the DY ratio 



of Eq. ( 64 ) . The ratio could be predicted from the 



known nucleon PDFs and the pion PDFs that have 
been extracted from this experiment (assuming no 
nucleon CSV) . The nucleon CSV distributions can 
then be extracted by comparing the predicted DY 
ratio with the observed value. Since the DY ra- 



tio of Eq. ( 64 ) results from subtracting two large 



and approximately equal terms, it is necessary to 
determine the relative DY cross sections to a few 
percent in order for this ratio to be statistically 
meaningful. Note that one can also exploit the 
fact that the CSV contribution to the DY ratio 
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depends only on x while the sea-valence term de- 
pends upon both x and x v . If the CSV term is suf- 
ficiently large, the process of extracting the CSV 
distributions may have to be carried out in an it- 
erative fashion. 



2. Parity-Violating Asymmetry in Electron 
Scattering 

The observation of parity violation in the scat- 
tering of polarized electrons from the deuteron, 
carried out in 1978 by Prescott et al. |139| . played 
a major role in validating the Standard Model. Re- 
cent advances in the technology of parity- violating 
experiments provide the possibility of repeating 
such experiments with an increase in precision 
of better than an order of magnitude |11U 1140] . 
These new experiments would allow a new pre- 
cision measurement of the Weinberg angle, they 
could probe physics beyond the Standard Model 
at the multi-TeV scale, and could provide tight 
constraints on nucleon parton distribution func- 
tions at large Bjorken x. In particular, the ratio 
d{x)/u{x) at very large x is not well determined 
|142j . As we will show, parity- violating electron 



scattering also has the possibility of observing par- 
ton charge symmetry violation at large Bjorken x. 

The parity violating (PV) asymmetry A PV for 
electron scattering on a nucleon can be written to 
lowest order in the 7 — Z interference in terms of 
the structure functions 



A PV {x,y) 

f(y) 
n(x) 

9% 



—G F Q 2 

1 - (i - vf 

[l + (l-y)2 



[g e A ri{x) + f(y)g e v r 2 (x)} 
E' 



Ff(x) _ 2E q e q g q v q + (x) 



Ff{x) 



Y,q e l<l + { X ) 



2i7(*) E q e 2 q q + (x) 
-l+4sin 2 6> w , g e A = -l 



In Eq. ( 66 ) , we have dropped some small correc- 
tions to the quantity f(y) and we have assumed 
the Bjorken limit where the longitudinal cross sec- 
tion is negligible relative to the transverse cross 
section. The additional terms are included in work 
by Hobbs and Melnitchouk [141 



In the proposed 
JLab PV experiment, the incident electron energy 
E will be in the range 10-11 GeV and the outgoing 
E' will run from 2 to 4 GeV. The parton model 



expressions for the ratios of structure functions are 
given by the quantities r± and r 2 in Eq. ( 66 ) . 



If we confine our attention to the region of x 
above 0.3 then the contribution to Eq. (66 1 from 
sea quarks should be quite small. Assuming that 
electron-deuteron scattering is given by the im- 
pulse approximation (the sum of scattering from 
proton plus neutron), and also assuming parton 
charge symmetry the expression for PV e — D scat- 
tering can be written 



—G F Q 2 



[4 + m<4] 



(67) 



In Eq. ( 67 1 , for couplings at tree level we have 
Ql ~ 5 

4 = ^ 

3 5 



(2# 



9t) 



9a 



(68) 



In Eq. ( 68 ) , the quark vector couplings are given 
in Eq. ( 27 1 , and the quark axial couplings are 



9a 



9 a 



(69) 



In this region, and with these assumptions, the PV 
asymmetry for e — D scattering depends weakly 
on y (the second term in Eq. ( 67 1 is significantly 



smaller than the first term) and is independent of 
x and of quark PDFs. 

We can now include the lowest-order CSV con- 
tributio n to the parity- violating e — D asymmetry. 
In Eq. ( 67 1 , the terms af and 03 are modified to 



d(0) 

n 

rf(0) 



S (csv) a 
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- Sd(x) _ 
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u(x) 
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Su(x) 


- Sd(x) 


af) 


.-10- 


2(2^ 


-9 d A )\ 


u(x) 


+ d(x) 
(70) 



(66) We note that in Eq. (70 1, the largest contribution 



to the CSV effect in the parity-violating electron 
scattering asymmetry comes from the CSV contri- 
bution to the denominator, i.e. from the structure 
function F^ D (x) (this is the origin of the 3/10 
term). The CSV terms will produce a correction 
to the PV asymmetry which has a characteristic 
dependence on Bjorken x. 

Fig. [20] plots the change in the e — D PV asym- 
metry 8A e J^/A e J^ arising from CSV effects, cal- 
culated by Hobbs and Melnitchouk |141j . This is 
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FIG. 20: The contribution from partonic CSV effects 
to the PV asymmetry for e— D scattering with incident 
electron energy 10 GeV, from Ref. QH]. The CSV con- 
tribution is given by Eq. ( 70 1 . Curves are labeled by 



the value of k from the phenomenological fit of valence 
CSV distributions determined from the MRST group, 
Ref. [23] and Eq. pBl. The best fit k = -0.2 and the 



90% confidence limits are represented by k = +0.65 



and k — —0.8. Upper figure: 
figure: Q 2 = 10 GeV 2 . 



= 5 GeV . Lower 



obtained from Eq. ( 70 1, vs. Bjorken x, for electron 
incident energy 10 GeV. The first graph plots the 
ratio for Q 2 — 5 GeV 2 and the second graph is 
for Q 2 = 10 GeV 2 . The CSV PDFs are obtained 
from the phenomenological global fit to high en- 
ergy data from the MRST group [23] . The valence 
CSV distributions were parameterized using the 
form of Eq. ( 35 1 . The three dashed curves repre- 
sent different values for the overall parameter k. 
One curve shows the best-fit value k = —0.2. The 
outer curves represent the values k = —0.8 and 
K = +0.65; these two values denote the 90% con- 
fidence limit for the valence CSV allowed in the 
MRST global fit. 

Within the 90% confidence limit, the predicted 
CSV contribution to the PV asymmetry tends to 



increase with increasing Bjorken x. The magni- 
tude of the CSV contribution increases with in- 
creasing Q 2 ; for Q 2 = 10 GeV 2 , the CSV con- 
tributions allowed within the MRST 90% con- 
fidence level range between roughly -0.025 and 
+0.03 at a value x — 0.7. Thus if experiments 
could achieve a precision of about 1% in the asym- 
metry, it should be possible either to observe ef- 
fects of partonic CSV in this experiment, or alter- 
natively to put strong constraints on upper lim- 
its for partonic charge symmetry violating effects. 
Note that our results are model-dependent, as the 
MRST group chose the particular functional form 
given in Eq. (35) for their partonic CSV PDFs. In 
addition, for simplicity MRST neglected the Q 2 
dependence of the CSV parton distribution func- 
tions in their global fits to high energy data. 



There is an additional uncertainty in the parity- 
violating asymmetries. Since our predicted PV 
asymmetry becomes significant only at large x, we 
need to account for the fact that the d/u ratio 
in the proton is rather poorly known in this re- 
gion. This uncertainty was studied by Botje [142] 
who extracted quark PDFs from a QCD analysis of 
combined HERA and fixed-target data. A precise 
determination of d/u at large x comes from the 
NMC measurements of muon DIS on proton and 
deuteron targets [55] [Ml IE]- However, in both 
cases the limit on accuracy is not the data but 
the theoretical understanding of the EMC effect 
in deuterium. For example, the covariant treat- 
ment of Fermi motion and binding by Melnitchouk 
and Thomas [61] show that, contrary to the con- 
clusions in the original paper, the SLAC data was 
consistent with the pcrturbative QCD predictions 
for the u/d ratio as x — * 1. For x > 0.4 the er- 
rors on the QCD predictions grow fairly rapidly. 
This occurs because the electromagnetic coupling 
is weighted by the squared charge of the quark 
flavor. 



An independent measurement of d/u in the pro- 
ton can be obtained by measuring W production 
in p — p collisions. In these reactions a W + tends 
to be produced by annihilation of a u quark from 
the proton and a d from the antiproton, while a 
W~ is produced from a d quark in the proton and 
a u from the antiproton. Because the u quark car- 
ries a larger momentum fraction than the d, W + 
production will tend to be boosted in the proton 
direction while W~ will be boosted in the antipro- 
ton direction. One measures a forward backward 



2!) 



asymmetry A(yi), where 



duction of charged pions from valence quarks, by 



A(yi) 



da{e+) _ da(e~) 
dyi dyi 



dcr(e+) 



da(£- 
dyi 



(71) 



In Eq. (71 1, yi is the rapidity of the lepton arising 
from decay of the W, and da^^/dyi is the differ- 
ential cross section for charged lepton production. 
This is a convolution of the cross section for W 
production, with the relevant W —* Iv decay dis- 
tribution [T48iri44] . 

The forward-backward asymmetries have been 
measured by the CDF [H5] and DO [HO [TTfj 
groups at the Fermilab Tevatron for p — p colli- 
sions at y/s — 1.96 TeV. The asymmetries tend to 
be particularly sensitive to the slope of the d/u ra- 
tio. The higher the rapidity, the larger the x range 
for which d/u can be studied |147j . 

The d/u ratio could also be determined from 
large- a; parity-violating electron scattering on hy- 
drogen, since the PV amplitude preferentially cou- 
ples to the down quark. This information could in 
principle fix the d/u ratio in the proton and elimi- 
nate some of the uncertainty in PV DIS reactions 
on deuterium. 



3. Charged Pion Leptoproduction from Isoscalar 
Targets 

In the preceding section it was pointed out that 
DY processes for charged pions on nucleons can 
test CSV because the tt^ 1 contain different valence 
antiquarks. For this reason, semi-inclusive pion 
production, from lepton DIS on nuclear targets, 
could be a sensitive probe of CSV effects in nucleon 
valence parton distributions |131l 1148] . The cross 
section for this process is given by [149 



1 da%{x,z) N Nh {x,z) 



(7 7v (x) dz 



N 



Nh _ 



e?^ (x)D?(z) . (72) 



E 



The quantity N 



Nh 



in Eq. (|72j) is the yield of 

and 



hadron h per scattering from nucleon N 
D'Z(z) is the fragmentation function for a quark 
of flavor i into hadron h. D'^z) depends on the 
quark longitudinal momentum fraction z = Eh/ 'v, 
where Eh and v are respectively the energy of the 
hadron and the virtual photon. 

For pion electroproduction on an isoscalar tar- 
get, charge symmetry relates the "favored" pro- 



Nf: v (x,z) = 4N?*-(x,z). 



(73) 



In Eq. (|73|, N?J V + (x, z) represents the yield of 



fav 



7r + per scattering from the deuteron, via the fa- 
vored mode of production (for ir + (tt~) produc- 
tion, the "favored" mode of charged pion produc- 
tion is from the target up (down) quarks). The 
HERMES collaboration at HERA (T5U] has mea- 
sured semi-inclusive pion production from hydro- 
gen and deuterium. 

Londergan et al. |131l 1148] showed that the ratio 
R A (x, z) is sensitive to parton CSV effects, where 
the ratio is defined by 



R A (x,z) 
R A {x,z) 
A(z) 

C A (z) 



N D * (x, z) _ N D7r+ {x,z) 



4A(z) 



A(z) 



N D * + (x,z) - N D ^(x,z) 
C A {z) [R os (x)+R sv (x,z)]; 



8(1 + A(z)) 
(l + 4A(z))(4 + A(z))' 



(74) 



Eq. (74 1 is evaluated at moderately large x, where 
the sea/valence ratio is small. It has been ex- 
panded to first order in the CSV nucleon terms 
and the sea quark distributions. The term R cs (x) 

,(x, z) is a sea- valence 
A CSV part 



is given in Eq. ( 64 ) , and R 



interference term given in Ref. [148J 
of the fragmentation function has been dropped as 
theoretical estimates suggest that this term should 
be very small |148] , 



The ratio R A of Eq. ( 74 1 has an overall factor 



that depends only on z; the normalization of the 
ratio is chosen to make this term close to one for 
moderate values of z. The remainder of this ratio 
contains two terms. The first term depends only 
on x, and is proportional to the nucleon valence 
CSV fraction; it is identical to the term R cs (x) de- 



fined in Eq. ( 64 1 , which could be measured in pion 



Drcll-Yan reactions. The final term in Eq. (74 1 
depends on both x and z; it is proportional to 
the sea quark contributions, and becomes progres- 
sively less important at large x. 

Fig. |2T] plots the ir-dependent contributions to 
the ratio R A (x, z) of Eq. (74) vs. x at fixed z = 0.4 
|131j . The solid (dot-dashed) curves show the non- 
strange (strange) sea contributions to the ratio. 
The strange quark contribution is negligible ex- 
cept at extremely low x. The long dash-dot and 
dotted curves show the contributions from quark 
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CSV contributions from the MRST global fit with 
k = —0.8 and +0.65, respectively; these repre- 
sent the 90% confidence limits for the MRST CSV 
PDFs. At this value of z, the coefficient C A (z) 
from Eq. (74 1 has a value very close to 1. The 
CSV terms are substantial only for large x > 0.4. 
The ratio requires precise experimental measure- 
ments of the x-dependence of R A (x, z) for fixed z. 
Note that this depends critically on the validity of 
the factorization hypothesis for the semi-inclusive 
yields, as given by Eq. (72 1. 



q 0.2 




21: 



Contributions of various terms to the ratio 

0.4. 



FIG 

R A (x, z) denned in Eq. |74l vs. x at fixed z 



Solid (dot-dashed) curves: nonstrange and strange sea 
quark contributions. Long dash-dot (dotted) curves: 
CSV contributions from Eq. (741, for k — —0.8 and 
k, = +0.65, respectively. Curves are calculated for 
Q 2 = 2.5 GeV 2 . 



For x > 0.4, the contributions from charge 
symmetry violating PDFs are substantial, and 
they rapidly become the dominant contribution 
at larger x. Thus, at the levels determined by 
the MRST global fit, it would appear that pre- 
cise measurements of charged pion production in 
semi-inclusive D1S electroproduction reactions on 
deuterium have the possibility of observing these 
isospin-violating effects, or they would be able to 
lower the current allowed limits on partonic CSV 
effects. Theoretically it would be possible to ob- 
serve such effects in measurements of e + D — > 
tt^ 1 + X at Jefferson Laboratory. However, the va- 
lidity of Eq. ( 74 ) requires that factorization (the 
fragmentation function for quarks into pions) be 
valid to within a few percent, ft would be neces- 
sary to demonstrate that factorization is obeyed 
to a very high degree, at energies available at Jef- 
ferson Lab. 

The validity of factorization has been checked 
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FIG. 22: Ratios of proton to deuteron semi-inclusive 
charged pion electroproduction cross sections on pro- 
ton and deuteron at fixed x = 0.32 as a function of z, 
from Ref. |151| . Solid (open) symbols reflect data after 
(before) subtraction of coherent p production events. 
Symbols are offset slightly in z for clarity. Top curve: 
sum of 7r + and ir~ cross sections (first equation in 
Eq. (76l); bottom curve: difference of cross sections 
4 76b ) . The hatched area in 



(second equation in Eq. 
the bottom curve indicates systematic uncertainties. 
The shaded bands represent a variety of calculations 
in both leading and next-to-leading order QCD from 
CTEQ [152] and GRV [85]. 



for semi-inclusive deep inelastic reactions at Jef- 
ferson Lab energies. Navasardyan et al. |151| mea- 
sured charged pion electroproduction from p and 
D with 5.5 GeV electrons. They were compared 
with a factorization hypothesis 



da ~ e 2 q q(x, Q 2 ) D q ^(z, Q 2 )e 



1 + Acos4> + bcos(2(j>) 
2^ 



bp* 



(75) 



With factorization, the semi-inclusive cross section 
appears as the product of a parton distribution 
function depending on x but not z, times a frag- 
mentation function for a quark to a pion that de- 
pends on z but not x. Assuming factorization, one 
can derive expressions for ratios of the pion elec- 
troproduction cross sections on protons and deu- 
terium, 



4u + (x) + d + (x) 



gp(7T + ) + Cr p (7T~) 
0-b(7T+) + 01,(71-) 

0>(tT + ) - 0p(7r~) 

(t d (tt + ) - cr D (7i—) 3(u v (x) + d y {x)) 



b{u+{x) 

4:U V (X) - 



-d+(x)) 
d v (x) 
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(76) 

Fig. [22] shows ratios of pion electroproduction 
cross sections on p and D at fixed x = 0.32 vs. z. 
From Eq. ( 76 1 , these linear combinations should 
be independent of z. Now, a number of assump- 



tions have gone into Eq. ( 76 1 ; in addition to factor- 



ization, this relation assumes parton charge sym- 
metry, neglects contributions from heavy quarks 
and also contributions from any pj- dependence 
of parton distributions. Nevertheless, to within 
about ten percent the ratios show very little de- 
pendence on z for z < 0.7. The deviation from 
these curves for z > 0.7 results from the N — ► A 
transition region. Furthermore, the shaded bands 
show the ratio that is expected from phenomeno- 
logical parton distributions from the CTEQ col- 
laboration |152j and from GRV [55]. The experi- 
mental results show that factorization is a rather 
good approximation at Jefferson Laboratory ener- 
gies. 

Despite the perhaps surprisingly good agree- 
ment with the factorization hypothesis in this en- 
ergy region, nevertheless factorization is not suffi- 
ciently accurate to carry out tests of charge sym- 
metry violation in pion electroproduction reac- 
tions, at current Jefferson Lab energies. The re- 
lations in Eq. (74 1 would be more reliable at a 



possible future electron-ion collider, where factor- 
ization should be assured to a high degree. As 
Eq. (74 1 was derived in lowest order QCD, it is 



necessary to check whether the results remain es- 
sentially unchanged in NLO. 



4- Test of Weak Current Relation 



F. 



W+Nq 



Or) = F. 



12 



jv 
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From Eq. (28 1, at high energies the F 2 struc- 



ture functions for charge-changing neutrino and 
antineutrino interactions on an isoscalar target are 
equal except for contributions from valence quark 
CSV, plus strange and charm quark terms, i.e. 



FY +No {x 



f? JV0 (^Q 2 ) 



2 ) _ J? W N 0{ 

x \Sdy(x) — Su v (x) + 2 (s~(x) — c~(x)j] 



(77) 

These cross sections might be measured at vari- 
ous experimental facilities. At a high energy elec- 
tron collider, weak interaction processes such as 
e~ p — > v e X are no longer completely negligible 
with respect to the electromagnetic process e~p — > 



reactions were made at HERA |153j . where pre- 
cise structure functions and parton distributions 
were measured for momentum transfers Q 2 > 100 
GeV 2 . Tests of parton charge symmetry would 
require collisions of electrons with an isospin-zero 
nucleus such as the deuteron. Then by compar- 
ing charge-changing weak interactions induced by 
electrons and positrons the quantity in Eq. (77 1 
could be measured. This might be feasible at a 
future electron-ion collider. 

Alternatively with very high energy neutrinos, 
the ratio of Eq. (77 1 could be measured by com- 



paring W boson production on an isoscalar target 
induced by neutrinos and antineutrinos. Theoreti- 
cal estimates of this process were made by Londer- 
gan, Braendler and Thomas |154j . One constructs 
the ratio 



R w {x) 



2 IF. 



W + D 



(x) 



pW~D 



{*)) 



FW +D {x)+F™- D (x) 
Sd v (x) — Su v (x) + 2(s~(x) — c~(x)) 



Rcsv(x) + R s (x) 



(78) 



At sufficiently high energies, the only quantities 
contributing to the ratio R\y are "valence" strange 
and charm distributions, or valence quark CSV 
terms. The ratio could also be checked for any 
isoscalar nuclear target, replacing the nucleon par- 
ton distributions by their nuclear counterparts. 

The top figure in Fig. [23] shows the theoret- 
ical CSV contribution, Rcsv(x) from Eq. (78 1. 
The dashed curve is calculated for Q 2 = 100 
GeV 2 , the dot-dashed curve for Q 2 = 400 GeV 2 , 
and the dash-triple dot curve for Q 2 = 10, 000 
GeV 2 . The quantity Rcsv(x) is predicted to be 
greater than 0.02 provided x > 0.4, using CSV 
estimates of Rodionov et al. Ref . |74j . All theoreti- 
cal calculations predict that in the valence region, 
5d v (x) is positive and Su v (x) negative, so their 
effects should add, producing several percent ef- 
fects at the largest values of x. The term R s of 
Eq. (78 1, proportional to the difference between 



strange quark and antiquark distributions (we ne- 
glect possible contributions from charm quarks), 
is shown in the bottom figure in Fig. |23| 



e X. Charged-current cross sections in e 
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As we mentioned in Sec. II B and expanded upon 
in Sec. |III C 2[ there are now experimental mea- 
surements from which one can extract the strange 
quark asymmetry s~ (x) — s(x)—s(x). These come 
from production of opposite-sign dimuon pairs in 
reactions initiated by v or v beams, from the 
CCFR and NuTeV groups [SUES]. The first mo- 
ment (s~ (x)) must be zero since there is no net 
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FIG. 23: Top: theoretical estimates of CSV contribu- 
tion Rcsv to the ratio Rw{x) of Eq. (78 1 vs. x, for 
various values of Q 2 . b) Bottom: theoretical estimates 
of s- s difference, R B (x) of Eq. (78}. From Ref. [154] . 



strangeness in the nucleon. This means that if 
there is a nonzero strange quark asymmetry it 
must have at least one node in x. The latest anal- 
ysis of these results by Mason [53] gives a positive 
value for the second moment S~ = (xs~(x)); this 
is in agreement with analyses of these experiments 
by the CTEQ group [64]. 

Note that these results have the opposite sign 
for the strange quark asymmetry from the cal- 
culations of Braendler et al. shown in Fig. |23) 
Those results were calculated in the framework 
of 'meson-cloud' models [55J ES EH SU EH 113 
extended to include strange quarks. In these 
models the nucleon fluctuates to a configuration 
N — > K + Y, where Y represents a A or X baryon. 
The s is associated with the virtual kaon produc- 
tion, while the s quark resides with the residual 
strange baryon [55] . Fig 24 shows the quantity 
s(x) — s(x) calculated using the model of Mel- 
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FIG. 24: The strange quark asymmetry s(x) — s(x) cal- 
culated in the meson-cloud model of Melnitchouk and 
Malheiro, Ref. |42j . The curves correspond to different 
values for the NKY form factor. The shaded region 
is an eetimate of the uncertainty in the calculations. 



nitchouk and Malheiro [42 . The curves show val- 
ues of s — s calculated using various values for 
the N — » KY form factor, and the shading rep- 
resents the uncertainty in the calculations. In the 
Melnitchouk-Malheiro calculation the s(x) — s(x) 
difference also has the opposite sign from the ex- 
perimental determination, although it should be 
noted that to within one standard deviation the 
experimental result is consistent with zero. For 
both the meson-cloud and experimental values for 
s(x) — s(x), the magnitude of the strange quark 
contribution to the quantity Rw in Eq- ( 78 1 is 



comparable to the contribution arising from CSV, 
although the ^-dependence of the two contribu- 
tions is quite different. 



IV. CHARGE SYMMETRY VIOLATION 
FOR SEA QUARKS 

For valence quark charge symmetry, several the- 
oretical models give quantitatively similar predic- 
tions for charge symmetry violating PDFs. Esti- 
mates of valence quark CSV by Sather [76] and 
Rodionov et al. [71] are in rather good agreement, 
and both the magnitude and shape of those va- 
lence CSV PDFs agree quite well with the best 
phenomenological global fit from the MRST group 
[2"5] , The situation is much different for charge 
symmetry in the sea quark sector. It is consid- 
erably more difficult to construct reliable theo- 
retical models to estimate sea quark CSV effects, 
and until recently the phenomenological situation 
was less certain. As we have seen, one problem is 
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that the tests of charge symmetry generally com- 
bine effects from heavy quarks in addition to CSV 
terms. At sufficiently large x the contributions 
from heavy quarks should become extremely small 
relative to valence quark CSV effects. However, 
at small x strange quark contributions are signifi- 
cant. Unless these contributions are known quite 
precisely it is difficult to determine the upper lim- 
its on parton CSV in the sea. 



A. Estimates of Sea Quark CSV 

The MRST group also searched for the presence 
of charge symmetry violation in the sea quark sec- 
tor [23]. They chose a specific functional form for 
sea quark CSV, dependent on a single parameter 
5, 



u n {x) = d?{x) 1 + 5 ■ 



d n {x) = u p (x) 



1-5 



(79) 



With the form chosen by MRST, the net momen- 
tum carried by antiquarks in the neutron and pro- 
ton are approximately equal; although this quan- 
tity is not conserved in QCD evolution, the change 
in momentum carried by antiquarks in the neutron 
was found to be very small in the kinematic region 
of interest. Using Eq. (79 1 to represent sea quark 



CSV effects, the MRST group performed a global 
fit to a wide array of high-energy data, where the 
coefficient 5 was varied to obtain the best fit. [The 
MRST group refer to the sea quark CSV param- 
eter as 5; however we use the notation 5 to avoid 
confusion with our definitions of partonic charge 
symmetry] 

Somewhat surprisingly, evidence for sea quark 
CSV in the MRST global fit was substantially 
stronger than for valence quark CSV. The x 2 t hey 



25 



obtain is plotted vs. the parameter 5 in Fig 
The best fit was obtained for 5 = 0.08, i.e., an 8% 
violation of charge symmetry in the nucleon sea. 
The x 2 corresponding to this value is substantially 
better than with no charge symmetry violation, 
primarily because of the improvement in the fit to 
the NMC n — D DIS data [HS [MJ E] when u n 
is increased. The fit to the E605 Drell-Yan data 
[155] was also substantially improved by the sea 
quark CSV term. 

Note that the MRST parameterization does not 
satisfy the "weak form" of charge symmetry as 
described in Eq. (20). The weak form of charge 




FIG. 25: [color online] The X 2 obtained by MRST, 
Ref. [23] for a global fit to high energy data of parton 
distribution functions including se a qu ark CSV with 



the functional form defined in Eq. (791. x 2 i s plotted 



vs. the free parameter 5 (which MRST label 5). 



symmetry would require that the first moments of 
sea quark CSV should vanish, i.e. 



(5u) = (5d) = 



(80) 



In fact, using the MRST parameterization of 
Eq. fl79l! for sea quark CSV and either the MRST 



or CTEQ phenomenological sea quark parton dis- 
tribution functions, then both of the first moments 
in Eq. (80 1 would be infinite. 



Bcncsh and Londergan used quark models to es- 
timate the magnitude of sea quark CSVj78j. They 
included sea quarks in quark model wavefunctions 
and attempted to calculate the sign and magni- 
tude of sea quark CSV in such models. Reasonably 
model-independent estimates have been made for 
valence quark CSV, but calculations of sea quark 
CSV require additional assumptions, and such cal- 
culations are likely to have substantial model de- 
pendence. Bcnesh and Londergan predicted very 
small CSV effects for antiquarks. They estimated 
that the fractional amount of CSV in the sea, Sq/q, 
should be at least an order of magnitude smaller 
than the corresponding fractional CSV effects for 
valence quarks. In quark model calculations, qual- 
itative arguments would suggest that sea quark 
CSV effects should be small. The relative mag- 
nitude of CSV effects will be given approximately 

by 



5j 



SM 

W) 



(81) 



where (M) is the energy of the lowest contributing 
intermediate states, and 5M is the mass difference 
for intermediate states related by charge symme- 
try. For antiquarks, the lowest energy states are 
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four-quark states, whose energy is roughly twice 
the energy of the lowest diquark states that con- 
tribute for valence quarks. The mass difference 
between charge symmetric four quark states is ap- 
proximately 6M w M n — M p = 1.3 MeV, or three 
times smaller than the mass difference for minority 
valence quarks. This naive estimate suggests that 
sea quark CSV effects should be roughly an or- 
der of magnitude smaller than for valence quarks. 
Cao and Signal carried out meson-cloud calcula- 
tions of sea quark CSV [52] ; they suggest that the 
bag model calculations of Benesh and Londergan 
neglect higher-order contributions that might be 
substantial. 

In their quark model estimates of sea quark 
CSV, Benesh and Londergan found that the sea 
quark CSV contributions Su(x) and Sd(x) tended 
to be roughly the same magnitude and to have op- 
posite sign; this is similar to the situation with the 
valence quark CSV. However, from Eq. (79 1 the 



phcnomenological sea quark CSV form assumed 
by the MRST group obeys the relation 



Su(x) + Sd(x) — — 5[d(x) — u(x)] 



(82) 



The quantity d(x) — u{x) has been measured by 



the E866 group |156l 1157] . who compared Drcll- 
Yan pD and pp experiments. Their measurements 
are in agreement with measurements of the same 
quantity at HERMES [T5S]. At HERMES this 
quantity was extracted from semi-inclusive DIS ex- 
periments of charged pion production from e — p 
and e — D reactions. We will discuss the measure- 



ments of this quantity in Sec. IV C 1 



B. Limits on Sea Quark CSV 

There are relatively few experimental limits on 
charge symmetry violation for sea quark parton 
distributions. As we mentioned in Sec. |IIIC 1| one 
can obtain rather strong experimental constraints 
on sea quark PDFs from experimental measure- 
ments of the W charge asymmetry in pp reac- 
tions [94] . Such measurements put some limits 
on the magnitude of charge symmetry violating 
sea quark PDFs |93| , but they primarily rule out 
very large CSV for sea quarks. As we have men- 
tioned previously, most tests of parton charge sym- 
metry have contributions from heavy quark PDFs. 
Strong tests of parton CSV will require rather pre- 
cise knowledge particularly of strange quark par- 
ton distributions. Since sea quark parton distribu- 
tions increase quite rapidly at small x, it is diffi- 
cult to separate contributions from partonic CSV 



in this region from effects due to strange quarks. 
The results obtained from the global fits of the 
MRST group to high energy data [23], discussed 
in the preceding section, appeared to show a sev- 
eral percent CSV effect. 

Perhaps the most promising area to search for 
sea quark CSV is to look for contributions to DIS 
sum rules [551 I159j . The lowest-order sum rules 
often involve integrals of parton distributions over 
all x. As we discussed in Sec. II B valence quark 



normalization requires that the first moment of the 
valence quark CSV terms vanish. Consequently, 
when one integrates parton distributions over all 
x, the only remaining charge symmetry violating 
contribution will arise from sea quarks. As we will 
discuss in the following section, sea quark CSV 
terms will contribute to various sum rules, par- 
ticularly the Gottfried sum rule [71] , and possibly 
also the Adler sum rule 169 . 



1. W Production Asymmetry at a Hadron Collider 

Production of W bosons resulting from the scat- 
tering of protons on an isospin-zero target repre- 
sents an area where, in principle, one could test 
parton charge symmetry. On an isospin-zero tar- 
get, e.g., the deuteron, we are interested in semi- 
inclusive reactions of the type p + D — > W + + X 
and p+D — > W~+X. We can then define the sum 
of W + and W~ cross sections and the forward- 
backward asymmetry 



<r s (x F ) 
A(x F ) 



da(x F ) 

a s (x F ) - a s (-x F ) 
a s (x F ) + a s (-x F ) 



d<j(x f 



dx e 



(83) 



In Eq. ( 83 1 the Cabibbo-favored terms in a s are 
invariant under the transformation x F — > — x F for 
Feynman x F — x\ — x-i. In the forward-backward 
asymmetry A(x F ), the only terms that survive are 
charge symmetry violating terms plus heavy quark 
terms in the Cabibbo-unfavored sector. 

As was discussed in Sec. |III CTj when the CCFR 
group performed its initial analysis |171 ITS] of the 
"charge ratio" R c defined in Eq. ( [52] ) , at small 
x the charge ratio appeared to deviate from one, 
with the deviation growing with decreasing x. It 
was pointed out by Boros et al. [T5J [5D] that, if 
the analysis of this data was accurate, a likely ex- 
planation of this discrepancy was a surprisingly 
large violation of charge symmetry in the nucleon 
sea quark PDFs. This CSV term in the sea quark 
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distributions was sufficiently large that one would 
expect the forward-backward asymmetry A(x F ) of 
Eq. (83) to be as large as several percent [20] , 

However, as was discussed in Sec. |IIIC 1| the 
low- a; discrepancy disappeared upon re-analysis of 
the CCFR experiment [22 . Londergan et at. then 
analyzed the prospects for W forward-backward 
asymmetry |160j using the (much smaller) sea 
quark CSV obtained in the MRST global fit to 
high energy data [23]. The asymmetries obtained 
were extremely small, generally less than one per- 
cent. The sea quark CSV terms were substantially 
smaller than suggested by the original CCFR anal- 
ysis. In addition, the strange quark contributions, 
which were originally much smaller than the CSV 
terms, were now no longer negligible, and they 
tended to cancel the CSV contribution. As a result 
of these newer results, it was concluded that these 
W production asymmetries no longer represent a 
promising prospect to determine charge symmetry 
violation in quark PDFs. 



2. Limits on Charge Symmetry Violation for Gluon 
Distributions 



It is possible that gluon distributions are differ- 
ent for proton and neutron. In fact, if sea quark 
distributions are charge asymmetric, then this 
would lead to a small charge symmetry violation in 
gluon distributions since sea quark and gluon dis- 
tributions are coupled through the DGLAP evo- 
lution equations [2H] [27] [25]. It is also conceiv- 
able that some other mechanism might give rise to 
charge symmetry violation in gluon distributions. 

Piller and Thomas |161] pointed out that CSV in 
gluon distributions might be probed through mea- 
surements of heavy quarkonium production. For 
example, if one considers J/if) production arising 
from nucleon-nucleon collisions, then the differen- 
tial cross section can be written as the sum of two 
terms 



d 2 o-(cc) 
dx F dM 2 



E- 

i,3 S 



1 



4M 2 



Vgg9i(xb)9j(x t ) 



l [q i {x b )qj{x t ) + qi(x b )qj(xt)} 



(84) 



In Eq. (84 1, the first term represents the contri- 



bution from gluon-gluon fusion and the second 
term is from quark-antiquark annihilation leading 



to heavy quarkonium. The corresponding quanti- 
ties a gg and a qq represent the subprocess cross sec- 
tions for gluon-gluon fusion and quark-antiquark 
annihilation, respectively. 

If one focuses on J/tp production in proton- 
neutron collisions then charge symmetry violation 
in either the quark or gluon distributions will pro- 
duce a forward-backward asymmetry in the re- 
sulting cross sections. Hence such a forward- 
backward asymmetry in J/ip production would 
be sensitive to partonic CSV. If we define the 
forward-backward asymmetry in terms of Feyn- 
man x F — Xb — xt, we obtain 

_ da(J/ip) da(J/if)) 
Aa pn (x F ) = — — \ XF — \^ XF , (85) 

then it is straightforward to show that 
(x F ) ~ d gg g(x t )Sg(x b ) 

+ d qq [u(x t )5d(x b ) + u(x t )Sd(xb) 

+ d(x t )Su(xb) + d(x t )Su(xb)] 

- [x b <-> x t ] . (86) 



In Eq. ((86), 5g(x) = g p (x) - g n (x). Piller and 
Thomas showed that CSV contributions in the 
gluon distribution and in sea quark PDFs would 
be most important at small x F while contributions 
from valence quark CSV should dominate at large 
values of x F . 

The E866/NuSea group has recently measured 
T production in reactions arising from 800 GeV 
protons on hydrogen and deuterium targets |162j . 
At these energies, the dominant contribution to T 
production comes from gluon-gluon fusion. In this 
case one expects the resonance cross-section ratio 



(j(p4 

Mp- 



D 

-P 



T) 



i 



5g{x t ) 

2g{x t ) 



(87) 



Thus CSV in gluon distributions would be man- 
ifested by a deviation of the T production ratio 



in Eq. (|87| from one. In Fig. [26]the E866/NuSea 
collaboration plot the cross section ratio for T pro- 
duction in pp and pD reactions vs. target x. These 
are plotted as the solid circles in Fig. [26] Within 
statistics they find no measurable deviation from 
one. This is in contrast to significant deviations 
which were seen for the ratio of pD and pp Drell- 
Yan processes (these are shown as the open squares 
in Fig. 26 1. The asymmetry in DY processes arose 
from differences between u and d distributions in 



the proton, as will be discussed in Sect. IV C 1 



From the T measurement we can put upper limits 
of roughly ten percent on possible CSV effects in 
gluon distributions. 
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FIG. 26: [color online] The ratio a{p + D)/2a(p + p) 
as a function of target x, for T resonance production 
cross sections for p + D and p + p reactions from the 
E866/NuSea Collaboration, Ref. [162]. Solid circles 
are the T production cross sections. For comparison 
the open squares are the corresponding ratios for the 
E866 Drell-Yan cross sections from Ref. 11571. 



Charge Symmetry Contributions to DIS 
Sum Rules 



Sum rules can provide extremely useful informa- 
tion on a single moment of parton distributions. 
If we choose the first moment of quark PDFs over 
x, we can invoke the quark normalization condi- 
tions on both the valence quark distributions and 
on valence quark CSV, as expressed in Eqs (17 1 



and (18 I, respectively. Consequently, terms which 
contribute to lowest-order sum rules will depend 
on integers that represent valence quark normal- 
izations, plus the first moments of antiquark distri- 
butions and sea quark CSV. Sum rules that mea- 
sure the first moments of various structure func- 
tions are extremely useful in that if one chooses 
appropriate linear combinations of structure func- 
tions, contributions from heavy quark CSVs will 
cancel out. This removes a major source of uncer- 
tainty, since for sea quarks the contributions from 
sea quark CSV and from heavy quark distributions 
are generally difficult to separate. 

Two sum rules, the Adler sum rule [53] and 
Gross-Llewellyn Smith sum rule [70], can be di- 
rectly related to linear combinations of quark nor- 
malization integrals. We will discuss the Got- 
tfried sum rule (GSR) [7T] in the following sec- 
tion. Unlike the Adler or Gross-Llewellyn Smith 
sum rules, the "naive" Gottfried sum rule expec- 
tation Sq — 1/3 is obtained only if one assumes 
both charge symmetry for parton distributions, 
and equality of light sea quark distributions in 
the proton sea, u p (x) — d p (x), or more precisely 
the equality of the first moment of these light sea 
quark PDFs. The expectation that the light quark 



sea distributions should be equal is often referred 
to as SU (2) flavor symmetry. This is an unfor- 
tunate connotation since the u p (x) and d p (x) dis- 
tributions are not related by any underlying dy- 
namical symmetry. However, if all of the light sea 
quarks were generated through gluon radiation, 
then one would expect the sea quark distributions 
to be identical except for effects due to light quark 
mass differences and small electromagnetic effects. 

As we have discussed in |III| theoretical expec- 
tations for CSV effects in parton distributions are 
expected to be no larger than a few percent. In 
Sect. IIIC we showed that current experimental 
upper limits on CSV are of the order of several per- 
cent for x < 0.4, and larger than 10% for x > 0.4. 
It would therefore be useful to construct sum rules 
which could in principle distinguish between CSV 
effects and effects arising from sea quark flavor 
asymmetry. In this Section, we will review the 
current status of the Adler, Gross-Llewellyn Smith 
and Gottfried sum rules, with particular atten- 
tion to the contributions from partonic CSV ef- 
fects. Good reviews of DIS sum rules up to about 
1996 can be found in the Handbook of Perturba- 
tive QCD [55] and in the review by Hinchliffc and 
Kwiatkowski [159 . We reviewed sum rules and 
parton CSV in detail in our previous review arti- 
cle [T5]. There has been little change in the sta- 
tus of the experimental Adler and Gross-Llewellyn 
Smith sum rules since the publication of that re- 
view article. 



1. Gottfried Sum Rule: Sea Quark Flavor and 
Charge Symmetry 

In the past fifteen years we have obtained a great 
deal of quantitative information on sea quark fla- 
vor asymmetry (d v (x) ^ u p (x)) in the nuclcon. 
Information on the first moment of these distribu- 
tions can be extracted from measurements of the 
Gottfried Sum Rule (GSR) [71], obtained from the 
difference of F2 structure functions from charged 
lepton DIS on neutrons and protons. The Got- 
tfried sum rule is also known as the valence isospin 
sum rule [159 . Using Eq. (34), we obtain 



S G 



1 F?{x)-lF{x) 

dx ■ 

x 

1 2 f~ 
3~3i 
2 



dx [dP(x) - u p (x) 
dx [4<W(x) + 5u(x)] . 
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If the nucleon sea is charge symmetric, and the 
first moment of the proton antiquark distributions 
are equal, then we obtain the "naive" expecta- 
tion S G — 1/3. Earlier measurements of the GSR 
[II [301 GSJ DM EM] obtained results that ap- 
peared to be less than the naive expectation of 
1/3, but with significant error bars. The first re- 
ally precise GSR value was obtained by the New 
Muon Collaboration (NMC) [H ISUj. The NMC 
result [891 K EES], S G = 0.235 ± 0.026, was more 
than four standard deviations below 1/3. Assum- 
ing charge symmetry, this implies a substantial ex- 
cess <P > u p . This effect is much larger than can 
be accommodated by perturbative QCD. Next-to- 
leading order (NLO) and next to next to lead- 
ing order (NNLO) QCD calculations predict very 
small effects |166j . However, the NMC result could 
be due to a combination of charge symmetry and 
flavor symmetry violating effects. 

Note that all three of the sum rules which we 
discuss here - the Gottfried, Adler and Gross- 
Llewellyn Smith sum rules - involve dividing the 
F 2 or xF% structure functions by x. This empha- 
sizes the contributions from small x. In all of these 
sum rules we invoke quark normalization condi- 
tions as given in Eqs. (17 1 and (18). These nor- 



malization conditions hold only after integration 
over all x. In reality one measures the sum rule 
down to some smallest value x m i n . In that case 
one must estimate the contributions from the un- 
measured region < x < x m i n - 

The GSR provides information only on the first 
moment of proton sea quark differences. It was 
proposed to make a 'direct' measurement of sea 
quark flavor asymmetry by comparing Drell-Yan 
processes initiated by protons, on proton and 
deuteron targets. This was suggested first by Eric- 
son and Thomas |167j in the context of the pionic 
explanation of the EMC effect, and later by Ellis 
and Stirling |168j . In the Drell-Yan [DY] process 
|123j hadronic collisions produce opposite sign lep- 
ton pairs with large invariant mass. The charged 
leptons are formed from the decay of a virtual 
photon arising from annihilation of a quark (anti- 
quark) in the projectile with an antiquark (quark) 
of the same flavor from the target. 

The experiments measure the ratios of Drell- 
Yan (DY) cross sections for incident protons on 
deuteron and proton targets. Assuming the valid- 
ity of the impulse approximation, the ratio is given 
by 



R UY {x 1 ,x 2 ) 



pD 
DY 

2<T PP 



d(x 2 ) - Sd(x 2 ) 
u(x 2 ) 



(89) 



The last line of Eq. ( 89 ) follows in the limit of large 
Feynman xf — x\ — x 2 assuming that 



d{x) 
u(x) 



0. 



(90) 



where x\ and x 2 are respectively the longitudi- 
nal momentum fractions carried by the projectile 
(target) quarks or antiquarks. If charge symme- 
try holds then from Eq. ( 89 1 , in the limit of large 
xp, the ratio of pD to pp Drell-Yan cross sections 
would directly measure the ratio of the down anti- 
quark to up antiquark distributions in the proton, 
at a given value of x 2 . 

Experiment NA51 at CERN [169J measured 
Drell-Yan processes for 450 GeV protons on proton 
and deuteron targets, obtaining a ratio u p /d p = 
0.51 ± 0.04(stai) ± 0.05(sysi) for a single aver- 
aged point (x) = 0.18. The E866 group at Fermi- 
lab |156l 11571 1170j compared Drell-Yan processes 
for 800 GeV protons on liquid hydrogen and deu- 
terium targets. Fig. [27] shows results from E866. 
For values x 2 < 0.2 the ratio is greater than one, 
and appears to decrease at higher values of x 2 , 
perhaps becoming less than one at x 2 ~ 0.3. Both 
the NA51 and E866 experiments show a substan- 
tial excess d p > u p at small x. This measurement 
was confirmed by subsequent semi-inclusive DIS 
(SIDIS) measurements, comparing yields of posi- 
tive and negative pions from scattering of energetic 
electrons on proton and deuteron targets at HER- 
MES [158]. The E866 group obtained the first 
moment of the light sea quark difference 



(d- 



0.118 ±0.012 



(91) 



Comprehensive review articles on light sea quark 
asymmetries have been produced by Kumano [83 
and by Garvey and Peng |171j . 



However, in order to obtain the result of Eq. (91 1 



the E866 group assumed parton charge symmetry. 
Eq. ( 89 1 shows that the Drell-Yan ratios could 



also have contributions from sea quark CSV ef- 
fects. One cannot extract the magnitude of fla- 
vor symmetry violating effects without assuming 
sea quark charge symmetry, as was emphasized by 
Ma [3U 1172] . Ma claimed that in principle one 
could explain the entire effect on the Gottfried sum 
rule from parton charge symmetry violation, even 
if parton flavor symmetry was exact. However, 
this would require sea quark CSV effects an order 
of magnitude larger than those obtained in the 
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FIG. 27: The ratio d p {x)/u p {x) as a function of target 
x, obtained from Drell-Yan cross sections for pp and 
pD reactions. Solid circles are data from the E866 
Collaboration, Ref. 11561 1157] ; the open square is the 
NA51 point, [169] . Curves are results of calculations 
using various phenomenological PDFs. 



FIG. 28: The quantity d p (x) — u p (x) as a function of 
target x. Solid circles: data from the E866 DY cross 
sections scaled to fixed Q 2 = 54 GeV 2 , Ref. [T561IT57] . 
Open squares: results from SIDIS measurements at 



HERMES, [T35], corresponding to (Q 2 ) = 2.3 GeV . 
Curves are calculations using various phenomenologi- 
cal PDFs. 



MRST phenomenological fit to high energy data 
[23j . A more natural source of the d — u difference 
was predicted by Thomas |173j . This incorporates 
effects of the pion cloud of the nucleon; the pro- 
ton predominantly emits a 7r + , which contains a 
valence d qua rk, leading to an excess d > u. 

we plot the quantity d p (x) — u p (x). 



In Fig. 



2<s 



The solid circles are the E866 DY points from 
Ref. USSJIHT] scaled to fixed Q 2 = 54 GeV 2 . The 
open squares are the results from the SIDIS mea- 
surements at HERMES |158j , which correspond to 
an averaged value (Q 2 ) — 2.3 GeV 2 . The curves 
are calculations using various phenomenological 
PDFs from GRV98 [HS] (dotted curve), MRST 
[174] (dashed curve) and CTEQ5M pS] (solid 
curve) . 

We first review results obtained for the Got- 
tfried sum rule with the MRST2001 [175] global fit 
PDFs with no isospin violation. In this parameter- 
ization, the sea quark distributions are obtained at 
a starting scale Qq — 1 GeV 2 and values at higher 
Q 2 can be obtained through DGLAP evolution. 
The MRST2001 fit obtains 



S G = 1/3- 2/'3(d-u) = 0.266 



(92) 



The MRST result for the Gottfried sum rule is 
just over la above the NMC value. The E866 pp 



and pD DY data essentially determine the MRST 
value for the d—u asymmetry in the proton. These 
results are basically consistent with the NA51 DY 
point |169j and with the HERMES semi-inclusive 
DIS measurements of positive and negative pion 
production in ep and eD scattering |158j . 

If one includes the phenomenological sea quark 
CSV effects obtained by MRST [23], then the con- 
tribution to the Drell-Yan ratio in the limit of large 
xf will be 



R DY (xi,x 2 ) 



1 



1.08d(x 2 ) 
u(x 2 ) 



(93) 



Thus, including the MRST sea quark CSV term 
would decrease the extracted sea quark flavor 
asymmetry by roughly 8% (one needs to take care 
since Eq. ( 93 1 is true only in the limit of large 
Feynman x). In principle, one could insert the sea 
quark distributions from the 2003 MRST global 
fit including CSV, and calculate the contribution 
from sea quark CSV to the Gottfried sum rule. 
However, using the MRST functional form from 
Eq. (79 1 gives an infinite result for S G . This is 



due to the fact that MRST choose the sea quark 
CSV proportional to the sea quark PDFs, which 
have infinite first moment. The MRST group 
plans to carry out future global fits assuming a 
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modified functional form for sea quark CSV [176] ; 
this would produce finite CSV contributions to the 
Gottfried sum rule. 



From Eq. ( 89 1 it is apparent that in comparing 



pp and pD Drell-Yan cross sections, one has con- 
tributions from both sea quark flavor asymmetry, 
i.e., d(x) 7^ u(x), but also from parton charge sym- 
metry violation. Peng and Jansen |177j pointed 
out that one can obtain information on d/u from 
measurements of W or Z production in pp colli- 
sions, which have no CSV contributions. For ex- 
ample, if one measures the ratio of W + and W~ 
production then one obtains 



R(x F ) 

R{x F )(x F =0) 
R(%f)(x f »0) 



_ d<r/dx F (p + p -> W + ) 



da/dx F (p + p 
u(x) d(x) 
d(x) u{x) 
u(xi) d(x 2 ) 
d(xi) u(x 2 ) 



W~) 



(94) 



In Eq. (94), the final two equations are true in 



the limit where one neglects strange quark contri- 
butions. Peng and Jansen showed that these Re- 
production ratios were quite sensitive to different 
phcnomenological predictions for light sea quark 
distributions. 



2. Adler Sum Rule 

The Adler Sum Rule [69 is given by the integral 
of the F 2 structure functions for charged current 
v and V DIS on the proton. The Adler Sum Rule, 
Sa, is defined (in the limit Q 2 — > 00) as 



dx 



F 2 w ' p {x,Q 2 )-F 2 w ^{x,Q 2 ) 
2x 



W+p, 



dx [<(z) - d p v (x) (1 - \V td \ 2 ) - s~(x)] 



1 



(95) 



We obtain the result S A = 1 if we neglect the term 
\V t d\ 2 ~ 1 x 10~ 4 . The Adler sum rule thus re- 
quires subtracting the F% structure function for an- 
tineutrinos and neutrinos on protons, and dividing 
by x (this emphasizes the contribution from very 
small x). The Adler sum rule then follows from the 
normalization of the valence quark distributions. 
As a consequence of the algebra of SU(2) charges, 
the Adler sum rule has no QCD corrections. Since 
the Adler sum rule involves measurements only on 
the proton, it has no CSV corrections. 



Another name for the Adler sum rule is the 
isospin sum rule |159j . Note: different overall nor- 
malizations for the Adler sum rule appear in the 
literature. An alternative normalization is a factor 
of two larger than ours [321 1159] . Our normaliza- 
tion agrees with that used by the WA25 experi- 
mental group |178] 1179] . 

The best experimental data to date are from the 
WA25 experiment [T7H1[T75], who used the CERN- 
SPS wide band neutrino beams in the BEBC H 
and D bubble chambers. Note that the WA25 
measurements involve neutrino CC measurements 
on neutrons (e.<?.,deuterons) and protons, and not 
v and v on protons, as given in the definition of the 
Adler sum rule, Eq. ( [95] ) . This was done because 
v beams generally have much higher fluxes than 
V. The WA25 experiment substituted neutrinos 
on neutron targets using the relation 



)=Fr P (x,Q*) 



+ 2x [s (x) — 8u{x) — 5d(x)] 



(96) 



as was discussed in Sect[Tj Thus, the WA25 group 
does not measure the integral S A of Eq. (951, but 
instead measures a different quantity 



S A = 



dx 



F 2 w+ "(x,Q 2 )-F 2 w+p (x,Q 2 ) 



2x 



(i 



dx [v%(x) — d%{x) — Su(x) — Sd(x)] 

(97) 



S A — / dx [Su(x) + Sd(x)\ 
Jo 



From Eq. ( 97 1 we see that the difference between 



the Adler sum rule and the measurements from 
WA25 involves the first moment of contributions 
from sea quark CSV. If charge symmetry is exact, 
or if the "weak form" of charge symmetry holds 
(see Eq. (20 1), then the Adler sum rule would be 
identical to what was measured by WA25, i.e., 

S A = Sa- 



Fig. 29 shows the experimental situation regard- 
ing the Adler sum rule. The experimental points 
are from the WA25 experiment [1781 1179] . The 
experimental data are shown for several values of 
Q 2 . The average value is S A = 1.01±0.08 (stat)± 
0.18 (syst). However, the total vN cross sec- 
tion used by the WA25 group is smaller than the 
presently accepted value [Ml ITOg] . If the WA25 
value is readjusted to fit this total cross section 
their result becomes S A = 1.08 ± 0.08 (stat) ± 
0.18 (syst). 
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FIG. 29: Experimental results for the "Adler sum rule" 
from the WA25 group, Ref. JT78lfT79] . Note that what 
is measured is given by Eq. ( |97| l and not the Adler sum 
rule of Eq. |(95|. 



The results show no significant Q 2 dependence. 
The large errors arise from the factor l/x in the 
integral, Eq. (95), which gives a heavy weighting 



to the data at small x. The paucity of data in this 
region and the relatively large error bars there give 
a large uncertainty in the sum rule value. There 
are currently efforts underway to develop new gen- 
eration neutrino experiments, with substantially 
higher fluxes than were available in the past. If 
these efforts come to fruition, it might be possible 
to test the Adler sum rule with v and V beams on 
a proton target. 

We can use the phenomenological sea quark 
CSV amplitudes determined by MRST [23] to esti- 
mate the CSV contribution to the WA25 measure- 
ment. Assuming that the quark normalization in- 
tegral is indeed one, then inserting the MRST sea 
quark function of Eq. (82 1 into the Adler sum rule, 
one obtains 



S A — S 



dx \5u(x) + Sd(x)] 
l+5(d p -u p ) = 1.008 . 



(98) 



The MRST result for sea quark CSV implies a 
difference of less than 1% between the Adler sum 
rule and the quantity measured in the WA25 ex- 
periment. Note that this result could be strongly 
model dependent because of the functional form 
assumed by MRST. 

One could also in principle measure a sum rule 
using antineutrino beams on protons and deu- 



dx 



Fr P (x,Q 2 )-FW-"(x,Q*) 
2x 



= S A + I dx [Su(x) + Sd(x)] 
Jo 

= l-S(d p -u p ) = 0.992 . 



(99) 



The last line of Eq. ( 99 ) holds if we assume the 



MRST result for sea quark CSV. Note that if 
the "we ak form" of charge symmetry holds (see 
Eq. ([201, then S A = S A and S A = S A . We will dis- 
cuss this more in Sect. IIV C 41 in connection with 
a "charge symmetry sum rule." 

One could also imagine measuring a similar 
quantity on a nucleus, rather than the proton. If 
one compares the integral of the F2 structure func- 
tion for a nucleus with Z protons and N = A — Z 
neutrons, then one would expect 



f M A /M 

S$= dx 



^F^~ A (x,Q 2 ) 



2x 



pW + A 



2x 



Z - N 
A 



(100) 



In Eq. |100| the quantity M A is the nuclear mass 
and M the nucleon mass, the structure functions 
Fi are normalized per nucleon, and we have as- 
sumed the impulse approximation. Nuclear effects 
in neutrino DIS have been evaluated at length by 
Kulagin and Petti [5JD [571 [55] . They calculate nu- 
clear modifications to structure functions arising 
from three general sources. The first category of 
effects, incoherent scattering from bound nucleons, 
tends to affect structure functions mainly at large 
Bjorken x. These are evaluated by Kulagin and 
collaborators through Fermi motion and nuclear 
binding effects. It can be shown that Fermi mo- 
tion and binding effects give zero contribution to 
the nuclear Adler sum rule. The second types of 
corrections tend to affect the structure functions in 
the region x ~ 0.1. These can arise from off- shell 
effects or from nuclear modifications of the meson 
cloud. The third type of corrections arise from 
coherent nuclear shadowing or anti-shadowing ef- 
fects. These predominantly affect the structure 
functions in the region of low x < 0.1. 

Because of the isovector nature of the Adler sum 
rule, nuclear pion corrections give zero contribu- 
tion. Kulagin and Petti |58j require cancellation 
of the isovector off-shell and nuclear shadowing 
corrections to the Adler sum rule. This provides 
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constraints on these corrections. Kumano and col- 
laborators [531 EH [55] have also made systematic 
considerations of nuclear corrections to structure 
functions and nuclear parton distributions. 



3. Gross- Llewellyn Smith sum rule 

The Gross-Llewellyn Smith (GLS) Sum Rule 
[70] is derived from the F$ structure functions for 
neutrinos and antineutrinos. This is also called 
the baryon sum rule |159j . If we sum the F$ struc- 
ture functions for neutrinos and antineutrinos on 



0L3 Sum Rule: dCFK Data at Q e = 3 OeV 2 



an isoscalar target, then from Eq. (33) we obtain 

? W + N 



So 



1 dx 
o 2a; 



xF? ' ly °(x) + xF 3 w N °{x) 



u v (x) + d v (x) + s (x) 
(du v (x) + 5dv(x)) 



b(n f ) 



2 

a s (Q 2 ) 

7T 

a s {Q 2 ) 



dx 



a(n f ) 



a s (Q 2 



AHT 



(101) 



The result S GLS = 3 follows from the normaliza- 
tion of the quark valence distributions. An iden- 
tical prediction would be obtained using either a 
proton or neutron target in the sum rule. The 
Gross-Llewellyn Smith sum rule holds only in lead- 
ing twist approximation, and only to lowest order 
in the strong coupling constant a s . Our expression 
for the GLS sum rule thus includes a QCD correc- 
tion (the term in square brackets in Eq. (|101[)), 



which was derived by Larin and Vermaseren 
and the quantity AHT represents a higher twist 
contribution |181j . 

As is the case for the Adler and Gottfried sum 
rules, the Gross-Llewellyn Smith sum rule requires 
that the structure function be divided by x in per- 
forming the integral. This gives a strong weighting 
to the small- a; region, such that as much as 90% 
of the sum rule comes from the region x < 0.1. 
Of the three sum rules we discuss in this review, 
the GLS sum rule is experimentally the best deter- 
mined. The most precise value has been obtained 



by the CCFR collaboration |182j , which measured 
neutrino and antineutrino cross sections on iron 
targets, using the quadrupole triplet beam (QTB) 
at Fermilab. A summary of experimental details 
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FIG. 30: Experimental results for Gross-Llewellyn 
Smith sum rule, Eq.fToT] from CCFR group, Ref. [TgS] . 



Squares: 1^3(1), sum of neutrinos plus antineutrinos, 
at Q 2 = 3 GeV 2 . Dashed curve: analytic fit to xFz. 
Diamonds: approximation to the integral S G ls{ x ) of 
Eq. 102 Solid line: fit to the integral Sgls(x). 



for precision measurements using high-energy neu- 
trino beams is given in the review article by Con- 
rad, Shaevitz and Bolton [38] . 



In Fig. 30 we show the CCFR measurements 
and the experimental values of xF 3 (x) (the sum 
of XF3 for neutrinos plus that for antineutrinos) 
vs. x. They obtain cross sections at several val- 
ues of x and Q 2 . The squares give the value 
of xFs(x) interpolated to an average momentum 
transfer Q 2 = 3 GeV 2 (this is the mean Q 2 for the 
lowest a>bin in the CCFR experiment, since the 
lowest x values contribute the greatest amount to 
the GLS sum rule). The dashed curve is the best 
fit to xF 3 of the form Ax b (l - x) c . The CCFR 
reported value for the sum rule at this Q 2 value 
is S GLS = 2.50 ± 0.018 {stat) ± 0.078 (syst). The 
GLS sum rule is therefore known to 3%. Because 
of the large contribution to the GLS sum rule from 
small x, one measures xF 3 at various values of x, 
and evaluates the integral 



1 dy 

2y 



yFr No (y) + vFr No (y) 



(102) 

The Gross-Llewellyn Smith sum rule is then ob- 
tained by estimating the limit 



Sols — lim^o S G ls{x) 



(103) 



The solid curve in Fig. 30 is S GLS (x). 

A theoretical value for the Gross-Llewellyn 
Smith sum rule requires evaluating the QCD cor- 
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Status of ^di^ Measurement 

FIG. 31: Experimental results for Gross-Llewellyn 
Smith sum rule, and their errors, for a series of ex- 
periments, in chronological order from top to bottom. 



rections. Calculations of the GLS sum rule in- 
clude next-to-leading order QCD corrections, us- 
ing a QCD scale parameter Aqcd = 213 ± 50 
MeV. With this scale parameter and NLO QCD 
corrections, one obtains a theoretical prediction 
S GLS = 2.63 ± 0.04 [183]. The theoretical pre- 
diction is two standard deviations above the ex- 
perimental value. In Fig. [31] we show the evolu- 
tion over time of the GLS sum rule value. The 
measurements shown are from the CDHS [184] . 
CHARM [T53], CCFRR [12], and WA25 [175] col- 
laborations. There are also two points from the 
CCFR measurements, the first using the Narrow 
Band Beam (NBB) neutrino data [1861 QUI and 
the second using the QTB data [182 from the Fer- 
milab Tevatron. 

The errors on the GLS sum rule are now at a 
level where the value of the strong coupling con- 
stant a s is a major source of error. The CCFR 
group now have data on xF$ over a wide enough 
range of Q 2 that, together with renormalized data 
from several other experiments, they may be able 
to evaluate the GLS sum rule without extrapola- 
tion for a large range of Q 2 values. This raises the 
hope that one can calculate the Gross-Llewellyn 
Smith sum rule as a function of Q 2 , and use the 
resulting Q 2 dependence of the sum rule to deter- 
mine a s (Q 2 ). The CCFR group has recently re- 
calculated both the GLS sum rule, and the strong 
coupling constant a s |188] . With data of this qual- 
ity over a large Q 2 range, it may be possible to 
use the Q 2 dependence to put constraints on the 
strong coupling constant. Additional information 
regarding this procedure can be found in the thesis 
by Seligman [IS] . 



x pW No (x^ which form the integrand for 
the GLS sum rule, are obtained by taking the 
difference between cross sections for neutrinos 
and antineutrino charged-current processes on 
isoscalar targets [IS] . In the limit of exact charge 
symmetry, the Fi structure functions exactly 
cancel in this subtraction, and only the F3 
structure functions survive. However, CSV effects 
make additional contributions to this integrand, 



i.e. using Eq. (28) 



3tt 
2G 2 M N E 

2 ; 



(da vN ° jdx - da" N ° /dx) 



xFy +N °(x,Q 2 ) + xF™- N °(x,Q 2 ) 



F. 



W + N 



{x,Q- z )-Fr No M Z ) 



u%{x) + d%(x) + 3s (x) — c (x) 



3 1 

-6uy(x) + -Sd v (x) 



(104) 



In addition to the light valence quark distribu- 
tions, Eq. (104 1 contains additional contributions 



from strange and CSV 'valence' quark distribu- 
tions. However, the CSV amplitudes have no ef- 
fect on the GLS sum rule value. Since quark va- 
lence distributions obey the normalization condi- 
tions of Eqs. (17 1 and (18 1, the contributions from 



valence strange and CSV terms must integrate to 
zero. Note, however, that the valence quark CSV 
effects contribute to the integral S GLS (x) at any 
finite value of x, and that the CSV effects vanish 
only upon integration over all x. 

Since the GLS sum rule is evaluated on nuclear 
targets (for the CCFR measurement, on iron), we 
need to consider nuclear modifications of structure 
functions and their potential effect on the GLS 
sum rule. Such an investigation has been carried 
out by Kulagin and Petti [55], using methods that 
were summarized in our discussion of the Adler 
sum rule (see Sect. IVC2I. As for the Adler sum 
rule, the Fermi motion and nuclear binding correc- 
tions have zero effect on the GLS sum rule. In the 
limit Q 2 — > 00, the off-shell and nuclear shadow- 
ing corrections also tend to cancel. In Fig. 32 we 
show the GLS sum rule as a function of Q 2 \Thc 
experimental points are the CCFR measurements 
on iron, vs. Q 2 . The curves are nuclear calcula- 
tions of Kulagin and Petti [58] for various nuclei. 
The dotted curve is for the nucleon and the solid 
curve for iron. This gives an idea of the magni- 
tude of nuclear corrections to the GLS sum rule, 
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FIG. 32: Nuclear corrections to the Gross-Llewellyn 
Smith sum rule, as a function of Q 2 , from Ref. |58| . 
Data points are from the CCFR measurements on iron, 
Ref. [182] . Dotted curve: nucleon results; dot-dashed 
curve: results for carbon; solid curve: results for iron; 
dashed curve: results for lead. 



and their Q 2 dependence. In addition there may 
be additional corrections from the pseudo-CSV ef- 
fects suggested by Cloet et al. |119j and discussed 
in Sect. ITlTDl 



4- A "Charge Symmetry Sum Rule" 



In Sect. IIV C II we showed that the Gottfried 
sum rule contains contributions from both charge 
symmetry violation and from sea quark flavor 
asymmetry. If sufficiently accurate experimen- 
tal data can be obtained, one could define sum 
rules which could differentiate between effects due 
to parton charge symmetry violation and those 
arising from differences in light sea PDFs, i.e., 
d p (x) ^ u p (x). Ma [35] defined a "charge symme- 
try" sum rule in terms of the F 2 structure func- 
tions for charged current neutrino and antineu- 
trino interactions on the neutron and proton, 



1 dx 



tpW + D 



F, 



W + p 



(x) 



7 W~p 



(x) 



{x)-F 2 w "(x) 



2 / dx [Su(x) + Sd(x)] 
'o 



(105) 



In Eq. ( |105| ), F 2 w+D (x) is the F 2 structure func- 
tion per nucleon for neutrino charged-currcnt DIS 
on the deuteron. From Eq. (105) we see that if ci- 



ther the strong form or weak form of charge sym- 
metry holds for the nucleon sea quark distribu- 
tions, then S C s w iU be zero. A deviation of this 
sum rule from zero would signal either a violation 
of parton charge symmetry, or a contribution from 
higher- twist terms. The higher twist terms would 
be expected to become progressively smaller with 
increasing Q 2 . Just as for the Adler sum rule, 
there are no QCD corrections to the charge sym- 
metry sum rule. 



The charge symmetry sum rule of Eq. (1051 is 



closely related to what has been measured as the 
Adler sum rule (see the discussion in Sect. IV C 2 1. 
We can easily see that 



dx 



F 2 w+ "(x,Q 2 )-F 2 w+p (x,Q 2 ) 



2x 



= S A 



~2~ 



(106) 



The WA25 group |178j measured cross sections 
from neutrinos on protons and deuterium, so their 
integral should give the Adler sum rule minus one- 
half the charge symmetry sum rule. However, as is 



shown in Fig. 29 errors in the WA25 measurement 
are of the order of 20%, so the charge symmetry 
sum rule at present is consistent with zero at the 
40% level, if we assume that the Adler sum rule is 
one. We can obtain an estimate of the charge sym- 
metry sum rule from the MRST global fit includ- 
ing sea quark CSV discussed in Sect. |IVB| Using 
the MRST form for sea quark CSV, Eq. (|98j) with 

best value 5 = 0.08 predicts S cs 0.016. The 

predicted quantity is extremely small; however, we 
stress that this result is based upon the (strongly 
model-dependent) functional form for sea quark 
CSV chosen by MRST (see Eq. (I79)). 



SUMMARY AND OUTLOOK 



We have reviewed the features of charge sym- 
metry, an approximate symmetry in particle and 
nuclear systems, as it relates to parton distribu- 
tions. First, we reviewed the relation between high 
energy cross sections and parton distributions, in 
terms of structure functions. Then we wrote the 
structure functions in terms of parton distribu- 
tion functions. Until recently, phenomenological 
parton distribution functions assumed the validity 
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of charge symmetry. This meant that PDFs for 
the neutron could be given in terms of those for 
the proton. If we relax this assumption, we must 
differentiate between neutron and proton PDFs. 
This requires the introduction of charge symmtry 
violating PDFs. 

In Sect. |II C| we expanded the structure func- 
tions in terms of parton distributions without 
making the assumption of charge symmetry. We 
then derived relations between these structure 
functions. Some relations that exist in the limit 
of exact parton charge symmetry must be modi- 
fied if we relax that assumption. 

In Sects. |III A| and |IIIB| we reviewed the phe- 
nomenological and theoretical situation regarding 
parton charge symmetry. One theoretical method 
for predicting partonic CSV contributions is to ex- 
amine the dependence of quark models for par- 
ton distributions on variations in quark and nu- 
cleon mass. For valence quarks, we showed that 
this leads to predictions for the magnitude and 
sign of the CSV terms. The quantities Su v (x) = 
u v p (x) — d v n (x) and Sd v (x) = d v p (x) — u v n (x) are 
found to be opposite in sign and roughly equal 
in magnitude. Since at large Bjorken x one has 
dv(a;) << u v (x), the fact that the valence CSV 
distributions are predicted to be roughly equal 
implies that the percent charge symmetry viola- 
tion for the "minority" valence quark distribution 
should be substantially greater than for the "ma- 
jority" valence quark distribution. 

There now exist phenomenological valence CSV 
PDFs from the MRST group [23]. They assumed a 
particular functional form with one overall free pa- 
rameter. That parameter was varied in a global fit 
of high energy experimental data. MRST obtained 
a very shallow minimum in fitting the high energy 
data. At the 90% confidence level the free pa- 
rameter k multiplying this phenomenological form 
could range between —0.8 and +0.65, with a best 
fit value k = —0.2. MRST chose a functional 
form such that Su v and dd v were required to be 
equal and opposite. This assured that the total 
momentum carried by valence quarks in the neu- 
tron and proton were equal (this quantity is rea- 
sonably well fixed by experiment). This choice by 
MRST agrees reasonably well with theoretical va- 
lence quark CSV PDFs obtained from quark mod- 
els. In fact the best fit of MRST is in surprisingly 
good agreement with theoretical valence CSV par- 
ton distributions from Sather |76j and Rodionov et 
al. UM- 

In Sect. HUB ll we discussed an additional mech- 
anism for charge symmetry violation. This oc- 



curs when a quark radiates a photon, in analogy 
with the well-known case where a quark radiates a 
gluon. If one incorporates these photon radiation 
terms into the QCD evolution equations, these 
give rise to a new type of charge-symmetry vio- 
lation. These "QED splitting" terms have been 
analyzed by two groups J2H [25] with qualitatively 
similar results. 



In Sect. Ill C we reviewed experimental limits on 



parton charge symmetry. Parton charge symmetry 
violation has never been directly observed. The 
strongest upper limits on parton CSV come from 
comparison of the F2 structure functions obtained 
from charged lepton DIS with those extracted from 
charged-current DIS arising from neutrinos and 
antineutrinos, with both taken on isoscalar tar- 
gets. We reviewed the upper limits that can be 
extracted by comparison of the NMC p — D reac- 
tions [H GEO GEE] with charged-current DIS for v 
and v on Fe, from the CCFR [T7] and NuTeV [22] 
experiments. Comparison of these experiments 
requires a number of corrections, from absolute 
cross sections normalizations to nuclear effects to 
shadowing corrections. However, the most thor- 
ough studies to date place upper limits of par- 
ton CSV at about the 6 — 10% level in the region 
0.03 < x < 0.4. 

The NuTeV group has obtained an independent 
measurement of the Weinberg angle by measur- 
ing charged-current and neutral-current DIS for 
v and v on an iron target [JJJ1 [3D]. They ob- 
tain a value for sin 2 9 W that differs by three stan- 
dard deviations from the value obtained at the Z 
pole. In Sect. |III C 2| we review this situation in 
considerable detail. We reviewed possible contri- 
butions from a number of "QCD effects" on the 
NuTeV measurement. In particular, we showed 
that valence parton CSV effects had the possibility 
to make substantial contributions to the NuTeV 
Weinberg angle measurement. At the 90% con- 
fidence level obtained in the MRST phenomeno- 
logical fit, partonic CSV could completely remove 
the NuTeV anomaly, or alternatively could make 
it twice as large. We showed that theoretical CSV 
effects predicted by quark models and obtained 
from QED splitting arc likely to remove approx- 
imately two-thirds of the NuTeV anomaly in the 
Weinberg angle. 

In Sect. IIIIDI we discussed a new nuclear reac- 
tion mechanism. This concerns the differential ef- 
fect of p exchange on protons and neutrons. Cloct 
et al. [119] have pointed out that this will produce 
effects that mimic those of CSV in a nucleus with 
N > Z. As a result we refer to this as pseudo 
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CSV. It was pointed out that this should produce 
a characteristic A dependence of the EMC effect; 
experiments have been proposed to look for evi- 
dence of this effect. It is predicted that this effect 
would account for up to two-thirds of the NuTeV 
discrepancy in the Weinberg angle. 

The current upper limits on partonic CSV for 
valence quarks are still reasonably large. Both the 
90% confidence level of the MRST phenomenologi- 
cal fit, and the best direct measurement from com- 
parison of the Fi structure functions for charged- 
lepton and charged-current DIS, give upper lim- 
its on valence quark CSV of the order of a few 
percent. Medium and high-energy facilities have 
now reached a precision where one could envision 
dedicated experiments that would either measure 
parton charge symmetry violation, or significantly 
reduce the upper limits on partonic CSV. 

In Sect. HIE we discussed four such experi- 



fects at small x. On rather general grounds one 
can argue that the magnitude of partonic CSV ef- 
fects should be given by 



ments. Since tests of parton charge symmetry 
require comparison of PDFs in the neutron and 
proton, all of these experiments require isoscalar 
targets. The first of these is a comparison of Drcll- 
Yan cross sections for n + and n~ projectiles on a 
target such as deuterium. Studies in the valence 
region for both pion and nucleon have the possibil- 
ity of revealing valence quark CSV effects. A sec- 
ond experiment would be parity-violating DIS in 
e — D reactions. Such an experiment is currently 
being planned for Jefferson Laboratory following 
the 12 GeV upgrade; it would have the potential 
to measure CSV effects if the PV asymmetry could 
be measured to roughly one percent. 

A third experiment that might test parton 
charge symmetry is a comparison of 7r + and n~ 
electroproduction from deuterium. Sufficiently 
precise experiments also have the possibility of re- 
vealing CSV effects. Identification of CSV effects 
requires that factorization be valid to a few per- 
cent, so such experiments might be most reliably 
carried out at a future electron-ion collider. One 
final experiment would be a comparison of W + 
and W~ production from neutrino and antineu- 
trino charged-current DIS on an isoscalar target. 
In all of these cases we have made estimates of the 
magnitude of CSV effects expected in these reac- 
tions. 

In Sect. [IV| we reviewed the situation regarding 
sea quark CSV effects. Theoretically, the situation 
for sea quark charge symmetry is not nearly as 
well founded as for valence quarks. In the valence 
region, sea quark contributions are quite small, 
whereas they are substantial at small x. It is very 
difficult to disentangle heavy quark and CSV ef- 
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In Eq. (107 1, the quantity dm would represent 



quark mass differences, in the range 1 and 5 MeV, 
and (M) denotes an effective mass of the system 
after removal of a quark. For valence quark CSV 
one would estimate (M) ~ 500 MeV, or a typ- 
ical diquark mass, while for sea quark CSV one 
would expect something like (M) ~ 1.3 GeV, a 
typical mass for a three quark-one antiquark state. 
From this rather general argument one would ex- 
pect that sea quark CSV effects should be signifi- 
cantly smaller than those for valence quark CSV. 

In Sect. |IV A"| we discussed the phenomenological 
sea quark CSV studies by the MRST group [53]. 
They assumed a functional form for sea quark CSV 
with an overall free parameter which was varied 
in a global fit to high energy data. Their best 
fit was obtained with a surprisingly large value, 
about 8%, for sea quark CSV. As we mentioned, 
at small values of x sea quarks, gluons and heavy 
quarks all contribute, making it difficult to isolate 
sea quark CSV effects. 

Since sea quark and gluon distributions are 
connected through QCD evolution equations, sea 
quark CSV should in principle lead to charge 
symmetry violation in gluon distributions. In 
Sect. |IV B 2| we reviewed possibilities for testing 
CSV in gluon distributions. A recent experiment 
measuring T production in pp and pD scattering 
could place upper limits on gluonic charge symme- 
try violation [162J. 

In Sect. |IV C| we reviewed one possibility to 
search for sea quark CSV, which is to look for con- 
tributions to DIS sum rules. Several of these sum 
rules involve the first moments of parton distribu- 
tions. The valence CSV distributions, and heavy 
quark "valence" distributions, must give zero first 
moment in order to respect valence quark normal- 
ization. In some of these sum rules the only terms 
whose first moment survives are sea quark CSV 
contributions. 

We showed that CSV contributions to the nu- 
cleon sea have no effect on the Gross-Llewellyn 
Smith or Adler sum rules, but in principle they 
affect the Gottfried sum rule. The Gottfried sum 
rule has contributions both from asymmetries in 
the light quark sea (the fact that d(x) ^ u(x)), 
and also from sea quark CSV. We also pointed out 
that the best existing experimental "test" of the 
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Adlcr sum rule actually measures a somewhat dif- 
ferent quantity, one that contains a non-zero con- 
tribution from sea quark CSV. The "Adler sum 
rule" measurement can thus be used to place an 
upper limit on parton sea quark CSV. Finally, we 
introduced a new sum rule, a "charge symmetry" 
sum rule which would be zero if either the "strong 
form" or the "weak form" of charge symmetry 
holds. A test of this charge symmetry sum rule 
would require measuring the structure functions 
for neutrino and antineutrino charged-current re- 
actions on protons and deuterium, with particular 
attention to the small- a; region. If one assumes the 
validity of the Adler sum rule then similar infor- 
mation could be obtained from measurements of 
either neutrinos or antineutrinos on an isoscalar 
target. 

In conclusion, in recent years much progress has 
been made in precision measurements of structure 
functions. From these one can extract parton dis- 
tribution functions which are now known to con- 
siderable accuracy. Recent experiments are now 
able to focus on specific questions such as the 
gluon distributions (both spin-independent and 
spin-dependent) and the flavor content of spin 
structure functions. We know that parton dis- 
tributions should have small charge symmetry- 
violating components. Recently, global fits of par- 
ton distributions have been carried out where one 
drops the assumption of charge symmetry. This 
gives indirect indications of the magnitude and 
shape of parton CSV. We have also discussed a 
series of experiments that could in principle reveal 



charge symmetry violation in parton distributions. 
Such experiments require great precision, coupled 
with an accurate knowledge of heavy quark dis- 
tributions. However, we are optimistic that such 
experiments can either lower the upper limits on 
parton CSV, or can find direct experimental evi- 
dence for parton charge symmetry violation. 
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